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1 Introduction 



In a sense this paper is a preliminary result of our recent attempt to analyse representations of the 
Sklyanin algebra. This is the algebra with four generators So, S a , a = 1,2, 3, subject to homogeneous 
quadratic relations 

[So, S a ]- = iJp^Sp, Sj] + , (1.1) 
[S a ,S f3 }-=i[S ,S 1 }+ (1.2) 



([A, B]± = AB ± BA, a triple of Greek indices a, (3, 7 in (1.1, |L2j) stands for any cyclic permutation of 
(1, 2, 3)). Structure constants of the algebra J a p have the form 



J a (3 — 



Jfi — Ja 



(1.3) 



where J a are arbitrary constants. Therefore, the relations (1.1 -Of) define a two-parametric family of 



quadratic algebras. These relations 
conditions under which operators 



1.1 



1.3) were introduced in the paper pj as the minimal set of 



L(u) 



3 

£ 

a=0 



W a (u)S a ®<T a , 



satisfy the equation 

R 23 (u - v)L 13 (u)L 12 (v) = L 12 (v)L 13 (u)R 23 (u - v). 

Here a a are Pauli matrices, Co is the unit matrix; W a (u) = W a (u\r], r), a = 0, . 
variable u with parameters rj and r: 

W a (u)= 9a+l{u) 



(1.4) 
(1.5) 

, , 3 are functions of the 



*a+lfa/2) (L6) 

(0 a (x) = 8 a (x\r) are standard Jacobi theta- functions with characteristics and the modular parameter r; 
for the definitions see Appendix to Section 6); 



R(u) 



a=0 



W a (u+^)a a 



(1.7) 



is the elliptic solution of the quantum Yang-Baxter equation 

R 23 {u - v)R 13 (u)R 12 (v) = R 12 (v)R 13 (u)R 23 (u - v) (1.8) 
that corresponds to the 8-vertex model. (The limit r — * yields the i?-matrix of the 6-vertex model.) 



I n (11.51) , (| 



we use the following standard notation. For any module M over the Sklyanin algebra 
eq. (1.4) defines an operator in the tensor product M®C 2 . The operator L 13 (u) (L 12 (u)) in the tensor 
product M ® C 2 <g> C 2 acts as L(u) on the first and the third spaces and as the identity operator on the 
second one (acts as L(u) on the first and the second spaces and as the identity operator on the third 
one). Similarly R 23 acts identically on M and coincides with the operator (1.7) on the last two spaces. 

Classification of discrete quantum systems solvable by the quantum inverse scattering method (see the 
reviews |p|j Q) is equivalent to solving eq. (1.5), where R(u) is a fixed solution to the Yang-Baxter 
equation fll.S|). 

More general elliptic solutions to (O]) were found in ||. Corresponding generalizations of the 
Sklyanin algebra were introduced in B, [Q. At present time only the simplest finite-dimensional repre- 
sentations of these generalized Sklyanin algebras are known. It would be very interesting to construct 
representations of these algebras in terms of difference operators similar to those found in || for the 
original Sklyanin algebra. 

As it was shown in the paper J8|, the operators S a , a = 0, ... ,3 admit representations in the form 
of second order difference operators acting in the space of meromorphic functions f(x) of one complex 
variable x. One of the series of such representations has the form 



(S a f)(x) = ^ 



Ja+l 



1+i (x - £rj)f(x + ??) - 9 a+1 (-x ~ £f])f(x - rj) 
01 (*) 



(1.9) 
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By a straightforward but tedious computation one can check that for any t, 77, t the operators (1.9) 
satisfy commutation relations (1.1-1. 3|) , the values of the structure constants being 



_ 9 a+1 (T,)8 a+1 (0) 

Therefore, the values of r and 77 parametrize the structure constants, while i is the parameter of the 
representation. Note that the original Sklyanin's parameter denoted by 77 in the paper S is equal to 
half of our 77 entering (JL6|), (Q, (jT|) and fll.lCQ . 



Putting /„ = /(rz.7j + xq), to the operators (1.9) we assign difference Schrodinger operators 



S a f n = A a n f n+1 +B a J n ^ (1.11) 

with quasiperiodic coefficients. The spectrum of a generic operator of this form in the space ^ 2 (Z) (square 
integrable sequences /„) has the structure of the Cantor set type. If 77 is a rational number, 77 = p/q, 



then operators (1.11) have q-periodic coefficients. In general (/-periodic difference Schrodinger operators 



have q unstable bands in the spectrum. 



In Sect. 5 of this paper we show that the spectral properties of the operator So given by eq. (1.9) 
are in this sense extremely unusual: 



Theorem 1.1 The operator So given by eq. \l.l\ ) for positive integer values of "spin" I and arbitrary 77 
has 2£ unstable bands in the spectrum. Its Bloch functions are parametrized by points of a hyperelliptic 
curve of genus 2£ defined by the equation 

2t+l 

y 2 = R(e)= l[(e 2 -ef). (1.12) 

z=l 

Bloch eighenfunctions ?/>(#, ±£j) of the operator So at the edges of bands span an invariant functional 
subspace for all operators S a . The corresponding At + 2-dimensional representation of the Sklyanin 
bra is a direct sum of two equivalent 21 + 1-dimensional representations of the Sklyanin 



Remark. In section 5 we show that there is a unique choice of signs for such that the space of irre- 
ducible representation of the Sklyanin algebra is spanned by Bloch eigenf unctions tp(x,£i). Unfortunately, 
at this stage we do not know at this stage any explicit constructive description of the corresponding split- 
ting of the edges in two parts. We conjecture that for real structure constants (when all e< are real) one 
has to choose all positive edges of bands: e$ > 0. 

This theorem indicates a connection between representations of the Sklyanin algebra and the theory 
of finite gap integration of solitonic equations. (The theory of finite difference Schrodinger operators B, 
fL0[ , ]ll[ , pa] was developed in the context of solving the Toda chain and difference KdV equations.) 
Furthermore, this theorem suggests that So is the proper difference analog of the classical Lame operator 

L = -^+£(£+l)p(x), (1.13) 

which can be obtained from So in the limit 77 — > 0. Finite gap properties of higher Lame operators (for 
arbitrary integer values of £) were established in p3[|. 



In Sect. 6 of this paper we suggest a relatively simple way to derive the realization (1.9) of the 
Sklyanin algebra by difference operators. This approach partially explains the origin of these operators. 
The basic tool is a key property of the elementary i?-matrix ([l]?]), which was used by Baxter in his 
solution of the eight- vertex model and called by him " pair-propagation thr ough a vertex" jL4j. A 
suitable generalization of this property for the arbitrary spin i-operator ( |l.4| ) leads to formulas (l.S). 



This derivation needs much less amount of computations than the direct substitution of the operators 



(1.9) into the commutation relations (1.1, 1.2). This method gives automatically the three representation 



series obtained by Sklyanin and an extra one which, presumably, was unknown. 

In the paper jl5| a remarkable connection between the motion of poles of the elliptic solutions of 
KdV equation (which are isospectral deformations of the higher Lame potentials) and the Calogero- 
Moser dynamical system was revealed. As it has been shown in Kg], fl7[|, this relation becomes an 
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isomorphism in the case of the elliptic solutions of the Kadomtsev-Petviashvili (KP) equation. The 
methods of finite gap integration of the KP equation were applyed to the Calogero-Moser system in the 
paper ]18| | , where the complete solution in terms of Riemann theta-functions was obtained. These results 
have been extended to spin generalizations of the Calogero-Moser system in the paper Jl!| . 

The main goal of this work is to extend this theory to elliptic solutions of the two-dimensional (2D) 
Toda chain and its non-abelian analogs. The equations of the 2D Toda chain have the form 

0+0_¥>„ = e *»-*»-i - evn+i-v-, d± = -^-. (1.14) 

Let us consider elliptic solutions with respect to the discrete variable n, specifically, 

¥>n = f(nrj + x , t+,i_) (1.15) 

such that 

c{x, t+,tJ) = cxp(Lp(x, t+,t-)- ip{x - 17, *+,*_)) (1.16) 

is an elliptic function of the variable x. In what follows we show that the function 
(exp ip) has the form 

exp<^(x,i + ,i_) = TT — X% + ^ , Xi=Xi(t+,t-), (1.17) 
fJi cr(x-Xi) 

(a(x\u>i, 0J2) is the standard Weierstrass a-function) and the dynamics of its poles x$ with respect to the 
time flows t +1 t_ coincides with the equations of motion of the Ruijsenaars-Schneider dynamical system: 
& 

ii = XiX s (V(xi - x s ) - V(x s - Xi)), (1.18) 

where 

V(x)=ax)-C(x + V), C(x) = Z tI' (1.19) 

a{x) 

This system is the relativistic analog of the Calogero-Moser model. Hamiltonians generating the com- 
muting i±-flows have the form 

«± = E e±w n ( " {x ' ~ VJ' - x r x ' ~ n) ) 112 < L20 > 

with canonical Poisson brackets: {pi,Xk} = 5ik- 

Our proof of this statement allows us to construct the action- angle variables for the system ( 1.18| ) 
and to solve it explicitly in terms of theta-functions. Being applyed to the non-abelian analog of the 2D 
Toda chain, this approach leads to spin generalization of the Ruijsenaars-Schneider model. 

This generalized model is a system of N particles on the line with coordinates Xi and internal degrees 
of freedom given by /-dimensional vectors = (a^ a ) and l-dimensional covectors bf = (bf) , a = 1, . . . , I. 
The equations of motion have the form 



5>+%)(&+a i )(nz« - xj) - V( Xj - Xl )), (1.21) 



J2aj{tfai)V(xi - Xj), (1.22) 



b+ = -^2b+(b+a j )V(x j -x i ). (1.23) 



The potential V(x) is given by the function (1.19) or its trigonometric or rational degenerations (V(x) — 
(cothx)^ 1 — (coth(a; + rj)) ,V(x) = x^ 1 — (x — ^)~ 1 , respectively). Hamiltonian formalism for this 
system needs special consideration and will not be discussed in this paper. 
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Let us count the number of non-trivial degrees of freedom. The original system has 2N + 2NI 
dynamical variables Xi, ii, cij, a , bf. The equations of motion (1.21-1.23) are symmetric under rescaling 



XiCii, hi 



Xi 



(1.24) 



The corresponding integrals of motion have the form ij = ii — (bfdi). Let us put them equal to zero: 

±i = (b+cn), (1-25) 



The reduced system is defined by N extra constraints J2 a bf = 1 (they destroy the symmetry ( 1.24 )). 
Therefore, the phase space of the reduced system has dimension 2NI. Moreover, the system has a further 
symmetry: 

a, W~ l Oi, b+ -> bfW, (1.26) 

where W is any matrix in GL(r, R) (independent of i) preserving the above condition on biS. This 
means that W must leave the vector v = (1, • • • , 1) invariant. Taking this symmetry into account, it is 
easy to see that dimension of the completely reduced phase space M is 



dim M. 



Nl 



(1.27) 



In the next three sections of this paper we derive explicit formulas for general solutions of the system 



(1.21-1.23) in terms of theta- functions. We would like to stress that these formulas are identical to those 
obtained for spin generalizations of the Calogero-Moser model in the paper |Q . At the same time the 
class of auxiliary spectral curves (in terms of which the theta- functions are constructed) is different. 
These curves can be described purely in terms of algebraic geometry. 

To each smooth algebraic curve L of genus N it corresponds a iV-dimensional complex torus J(T) 
(Jacobian of the curve). A pair of points s L defines a vector U in the Jacobian. Let us consider 
a class of curves having the following property: there exists a pair of points on the curve such that the 
complex linear subspace generated by the corresponding vector U is compact, i.e., it is an elliptic curve 
£q. This means that there exist two complex numbers 2u a , Im w^/wi > 0, such that 2u> a U belongs to 
the lattice of periods of holomorphic differentials on T. From pure algebraic-geometrical point of view the 
problem of the description of such curves is transcendental. It turns out that this problem has an explicit 
solution and algebraic equations that define such curves can be written as a characteristic equations for 
the Lax operator corresponding to the Ruijsenaars-Shneider system. Moreovere, it turns out that in 
general position £q intersects theta-divisor at N points a;, and if we move So in the direction that is 
defined by the vector V + {V~) tangent to T e J(T) at the point P + (-P - ), then the intersections of So 
with the theta-divisor move according to the Ruijsenaars-Shneider dynamics. An analogous description 
of spin generalisations of this system is very similar. The corresponding curves have two sets of points 
P i , i — such that in the linear subspace spanned by the vectors corresponding to each pair 

there exist a vector U with the same property as above. 

The following remark is in order. The geometric interpretation of integrable many body systems of 
Calogero-Moser-Sutherland type consists in the representation of the models as reductions of geodesic 
flows on symmetric spaces pl| . Equivalently, these models can be obtained from free dynamics in a larger 
phase space possessing a rich symmetry by means of the hamiltonian reduction ]22j . A generalization 
to infinite-dimensional phase spaces (cotangent bundles to current algebras and groups) was suggested 
in [p3[ , ]24| . The infinite-dimensional gauge symmetry allows one to make a reduction to finite degrees 
of freedom. Among systems having appeared this way, there are Ruijsenaars-Schneider-type models and 
the elliptic Calogero-Moser model. 

A further generalization of this approach should consist in considering dynamical systems on cotan- 
gent bundles to moduli spaces of stable holomorphic vector bundles on Riemann surfaces. Such systems 
were introduced by Hitchin in the paper |25[| , where their integrability was proved. An attempt to iden- 
tify the known many body integrable systems in terms of the abstract formalism developed by Hitchin 
was recently made in |2q] . To do this, it is necessary to consider vector bundles on algebraic curves with 
singular points. It turns out that the class of integrable systems corersponding to the Riemann sphere 
with marked points includes spin generalizations of the Calogero-Moser model as well as integrable 
Gaudin magnets [E7| (see also |28|). 
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However, Hitchin's approach, explaining the algebraic-geometrical origin of integrable systems, does 
not allow one to obtain explicit formulas for solutions of equations of motion. Furthermore, in general 
case any explicit form of the equations of motion is unknown. We hope that the method suggested for 
the first time in Jl8[ ) (for the elliptic Calogero-Moser system) and further developed in the present paper, 
could give an alternative approach to Hitchin's systems; may be less invariant but yielding more explicit 
formulas. We suspect that this method has not yet been used in its full strength. Conjecturelly, to each 
Hitchin's system one can assign a linear problem having solutions of a special form (called double-Bloch 
solutions in this paper). In terms of these solutions one might construct explicit formulas solving the 
initial system. 

Concluding the introduction, we remark that this paper can be devided into three parts which are 
relatively independent. The structure of the first one (Sects. 2 - 4) is very similar to that of the paper 
|^9| . Furthermore, in order to make our paper self-contained and to stress the universal character of the 
method suggested in |Q , we sometimes use the literal citation of the paper (l^] . At the same time we 
skip some technical detailcs common for both cases, trying to stress the specifics of difference equations. 
In the second part (Sect. 5) discrete analogs of Lame operators are introduced and studied. Finally, 
in the third part (Sect. 6) we give a simple derivation of difference operators representing the Sklyanin 
algebra and explain the origin of these operators. Actually, we expect a deeper connection between the 
three main topics of this paper; a discussion on this point is given in Sect. 7. 



2 Generating linear problem 



The equations of the non-abelian 2D Toda chain have the form 

d+dd-g^gn 1 ) = 9n9n-l - Qn+ign 1 - 



(2.1) 



These equations are equivalent to the compatibility condition for the overdetermined system of the linear 
problems 

d + ip n (t+,t-) = i/; n+1 (t + ,t-) +v n (t+,t-)ip n (t+,t-), (2.2) 
fl_^ n (t +J t_) = c n (i +J *_)Vn-i(*+,*-). (2-3) 

where 

Cn = g n g~-i, v n = (d+g^g' 1 (2.4) 



(g n is I x I matrix) . Similarly to the case of the Calogero-Moser model and its spin generalizations |L8| , 
p9| , the isomorphism between the system (1.21-1.23) and the pole dynamics of elliptic solutions to the 
non-abelian 2D Toda chain is based on the fact that the auxiliary linear problem with elliptic coefficients 
has infinite number of double-Bloch solutions. 

We call a meromorphic vector-function f(x) that satisfies the following monodromy properties: 

f{x + 2u a ) = B a f(x), a =1,2, (2.5) 

a double-Bloch function. Here uj a are two periods of an elliptic curve. The complex numbers B a are 
called Bloch multipliers. (In other words, / is a meromorphic section of a vector bundle over the elliptic 
curve.) Any double-Bloch function can be represented as a linear combination of elementary ones. 

Let us define a function <fr(x, z) by the formula 



z) = 



a{z + x + -q) 
a(z + rj)a(x) 



a(z - rj) 



-i x/2-q 



a{z + rj) 



(2.6) 



Using addition theorems for the Weierstrass tr-function, it is easy to check that this function satisfies 
the difference analog of the Lame equation: 



where 



$(x + r),z) + c(x)$(x —T),z) = E(z)<5>(x, z), 
a(x — rj)a(x + 2rf) 



c(x) 



a(x + r/)a(x) 



(2.7) 
(2.8) 



6 



Here z plays the role of spectral parameter and parametrizes eigenvalues E{z) of the difference Lame 
operator: 

E(z) = a( - 2 ^ (2.9) 
(a(z - rj)a(z + 77)) 1 / 2 

The Riemann surface Tq of the function E(z) is a two-fold covering of the initial elliptic curve Tq with 
periods 2u a , a = 1, 2. Its genus is equal to 2. 

Considered as a function of z, <&{x, z) is double-periodic: 

$(x,z + 2w Q ) = $(x,z). (2.10) 

For values of x such that a;/2^ is an integer, <E> is well defined meromorphic function on To. If x/2r) is 
a half-integer number, then $ becomes single-valued on IV For general values of x one can define a 
single- valued branch of <&(x, 2) by cutting the elliptic curve Tq between the points z = ±77. 

As a function of a; the function $(2;, z) is double-Bloch function, i.e. 

$(a; + 2w Q ,z) =T Q (z)$(x,z), (2.11) 

where Bloch multipliers are equal to 

T a (z)=eM^CM(z + v))(^^) ■ (2-12) 

In the fundamental domain of the lattice defined by 2ui a the function $(x, z) has a unique pole at the 
point x = 0: 

$(z,2r) = -+A + 0(ar), A = C(z + t?) + j- In ^ ~ T? j ■ (2-13) 
X 2r? cr(z + 77) 

That implies that any double-Bloch function f(x) with simple poles at points Xi in the fundamental 
domain and with Bloch multipliers B a such that at least one of them is not equal to 1 may be represented 
in the form: 

JV 

f(x) = J2 s M x - x i> z ) kX/n > (2-14) 

where z and a complex number k are related by 

B a =T a {z)k 2 ^/\ (2.15) 

(Any pair of Bloch multipliers may be represented in the form ( 2.15| ) with an appropriate choice of 
parameters z and k.) 

Indeed, let Xi, i = 1, . . . ,m, be poles of f(x) in the fundamental domain of the lattice with periods 
2wi,2o->2- Then there exist vectors Sj such that the function 

m 

F(x) = f(x) -J2 S M* ~ Xi,z)k x ^ 

i=l 

is holomorphic in x in the fundamental domain. It is a double-Bloch function with the same Bloch 
multipliers as the function /. Any non-trivial double-Bloch function with at least one of the Bloch 
multipliers that is not equal to 1 has at least one pole in the fundamental domain. Hence F = 0. 

The gauge transformation 

f(x) /» = f(x)e ax , (2.16) 

where a is an arbitrary constant does not change poles of any functions and transforms a double-Bloch 
function into another double-Bloch function. If B a are Bloch multipliers for /, then the Bloch multipliers 
for / are equal to 

B 1 = Bxe 2aUl , B 2 = B 2 e 2aul2 . (2.17) 



Two pairs of Bloch multipliers are said to be equivalent if they are connected by the relation (2.17) with 
some a. Note that for all equivalent pairs of Bloch multipliers the product 

B ui B -"i _ B (2.18) 
is a constant depending on the equivalence class only. 
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Theorem 2.1 The equations 

d t V(x, t) = V(x + r), t) + a,i(t)bf{t)V(x - xi(t))*(x, t), 



N 



i=l 



N 



d t V + (x, t) = V+(x- v , t) + V + (x, t) cii(t)b+{t)V(x ~ x t (t)) 



(2.19) 



(2.20) 



have N pairs of linearly independent double-Bloch solutions ^/^(x,t), tyfJXft) with simple poles at 
points Xi(t) and (xi(t) — rj), respectively, 



^ {s) (x + 2uj a ,t) = B a ^ (s) (x,t), V+ s) (x - 2w a ,t) =B a<s ^ (s) (x,t), 
with equivalent Bloch multipliers (i.e. such that the value 

BZ B 2? = B- 



(2.21) 



(2.22) 



does not depend on s) if and only if Xi(t) satisfy equations ( 1.2% ) and the vectors <2j, bf satisfy the 
constraints (1.25) and the system of equations 



di = a, J (b^a i )V(x l - Xj) - \ t 



(2.23) 
(2.24) 



where Xi — Aj(i) are scalar functions. 



Remark. The system fll.2l| ), fl2.23| ), fl2.24| ) i s "ga u ge eq ui valen t" to the system fll.2l| )-([l.23|). This 
means that if (xi,ai,bf) satisfy the equations (1.21), ( [2.23 ), (2.24), then Xi and the vector-functions 



(2.25) 



fij = diqi, b~\ = b i q i l , q l = cxp( / \i(t')dt') 
are solutions of the system ( |l.2l| - |h23 ). 



Theorem 2.2 If the equatio ns ( 2.1i ), (2.2C) have N linearly independent double-Bloch solutions with 
Bloch multipliers satisfying ( 2. 2i ) then they have infinite number of them. All these solutions can be 
represented in the form: 

N 

* = eft, k, z)9(x - Xl (t),z)k x /'\ (2.26) 



JY 



=J2sf(t,k,z)${-x + Xi(t)-r),z)k- x/7 >, 



(2.27) 



where Si is l-dimensional vector s$ = (sj )Q! ), sf is l-dimensional covector sf — (sf). The set of corre- 
sponding pairs (z, k) are parametrized by points of algebraic curve defined by the equation of the form 



N 



R(k,z) = k N + J2n(z)k 



N-i 



i=l 



Proof of Theorem 2.1. As it was mentioned above, ^i s \(x, t) (as any double-Bloch function) may be 
written in the form (2.26) with some values of the parameters z s ,k s . The condition (2.22) implies that 
these parameters can be choosen as follows: 

z s = z, s = l,...,N. 
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Let us substitute the function ^(x,t, z, k) of the form (2.26) with this particular value of z into eq. 
( [2.19 ). Since any function with such monodromy properties has at least one pole, it follows that the 
equation ( 2.19| ) is satisfied if and only if the right and left hand sides of this equality have the same 
singular parts at the points x = Xi and x = Xi — r\. 

Comparing the coefficients in front of (x — Xi)~ 2 in ( [2.19] ) gives the equalities 

x i s l = a,i(bfsi), (2.28) 

hence the vector Sj is proportional to ac 

Si t a(t, k, z) = a(t, k, z)ai ta (t), (2.29) 

and Qj, bi satisfy the constraints ( |l.25| ). Cancellation of the coefficients in front of (x — x% + r/)^ 1 gives 
the conditions 

— ksi + dibf Sj$>(xi — Xj —r],z) = 0. (2.30) 



Taking into account (1.25) and (2.29), we rewrite them as a matrix equation for the vector C = (q): 

(L(t, z) - M)C = , (2.31) 
where I is the unit matrix and the Lax matrix L(t, z) is defined as follows: 

Lij(t, z) = (bfaj)<Z>(xi - xj - r), z). (2.32) 
Finally, cancellation of the poles (x — xt)^ 1 gives the conditions 



k - [ ^2 a J b+V{x l - Xj ) + (A - t(rj))aibf s, - a< ^2(bfs j )^>(x l - Xj,z) = 0. 



(2.33) 



Using ( |2.29[) , we get the equations of motion ( [2. 23] ), where 

h{t) = - + (Cfa) - A)±i - ^(b+aiWxi - Xj ,z)- 



Rewriting (2.34) in the matrix form, 

(d t + M(t, z))C = , 

we find the second operator of the Lax pair: 

Mij(t, z) = (-Xi + (C(ri) - A)x i )5 ij - (1 - 5 ij )bja j ^(x l - Xj,z). 

Substituting in the same way the vector ^ + ( |2.27 ) into the equation ( 2.20| ), we get: 

sf(t,k,z) = cf(t,k,z)bf(t), 

C + {L{t,z) -kl) = 0, 
-d t C+ + C + M {+) (t,z) = 0, 



where C + = (cf), L is given by ( 2.32 ), and M^ + > has the form ( 2.36 ) with changing A,(t) to 



Kit) = + (C(tj) - A)±i - ^(bjaiMxj - x u z) 

Moreover, it follows that the covector bf satisfies the equations 

h t = - E *$<f>t<h)V{xi - Xi) + Xfb+ . 



(2.34) 

(2.35) 
(2.36) 

(2.37) 

(2.38) 
(2.39) 

(2.40) 
(2.41) 
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The assumption of the theorem imply that the equations ( [2.3l| ), (2.35) and (2.38), (2.39) have N linear 
independent solutions (corresponding to different values of k). The compatibility conditions of these 
equations have the form of the Lax equations 



L+[M,L]=0, L+[M (+ \L] = 0. 

In particular, they imply that Xi = Xf . 

The function ^(x, z) satisfies the following functional relations: 

$(x - T], z)<S>(y, z) - $(x, z)$(y - rj, z) = <&{x + y - rj, z)(V{-x) - V(-y)), 

&(x — r),z) — - r],z)(V(—x) + Civ) - A) - $>(—T],z)$(x,z) 



(2.42) 



(2.43) 
(2.44) 



(the constant A is defined in (2.13)). The first relation is equivalent to the three-term functional equation 
for the Weierstrass a- function. The second relation follows from the first one in the limit y — > 0. 



The identities ( 2.43 ) and ( 2.44 ) allow one to prove by direct computations the following lemma 
completing the proof: 



Lemma 2.1 For th e matrices L and M define d by ( 2. Si ) and ( 2.3t) resp ectively, where ai andu[ satisfy 
the equations ( 2. 25 ) and the constraints ( 1.21 ), the Lax equations ( 2.4-i ) hold if and only if x% (t) satisfy 
the equations (1.21). 



Proof of Theorem 2.2. As it was proved above, eqs. ( |2.19 ), ( |2.20|) hav e N linearly independent solutions 
if eqs. ( |2.23 ) and (2.24), constraints (3.£) and Lax equations (2.42) are fulfilled for some value of the 
spectral parameter z. But according to the statement of Lemma 2.1, the Lax equations are then fulfilled 
for all values of the spectral parameter z. Therefore, for each value of z we can define the double-Bloch 
solutions of eq. (2. IE) w ith t he hel p of formulas (2.26) and ( [2.29 ), where Cj are components of the 
common solution of eqs. ( [2.31 ) and ( [2.35 ). 

As it follows from (|2.3l|) , all addmissible pairs of the spectral parameters z and k satisfy the charac- 
teristic equation 

R(k, z) = det (kl - L(t, z)) = 0. 

At the begining of the next section we prove that this equation defines an algebraic curve T of finite 
genus. That will complete the proof of the theorem. ■ 



Remark 1. In the abelian case (I — 1) the operators of the Lax pair by a "gauge" transformation with 
the matrix Uij = aidij can be represented in the form 



LVj X) = Xi$(Xi - Xj - T), z), 



(2.45) 



M 



(7=1) 



(C(n) - A)±i - V(xi - x s )±s 5ij - (1 - 6ij)xi$(xi - Xj,z). 



(2.46) 



These formulas are equivalent to the Lax pair obtained in E{| , [Q ■ 

Remark 2. In the abelian case it is chough to require that only one of eqs. (pT9| ) (EE20l) has N linearly 



independent double-Bloch solutions with Bloch multipliers satisfying conditions (2.22) 



3 The direct problem 



As it follows from the Lax equation ( 2.42| ), the coefficients of the characteristic equation 

R(k,z) ee det (kI-L(t,z)) = (3.1) 



do not depend on time. Note that they are invariant with respect to the symmetries ( 1.24 ), ( |l.26| ) 
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Theorem 3.1 The coefficients ri(z) of the characteristic equation (3.1) 



N 



R(k,z) = k N + J2n(z)k N -* (3.2) 

i=l 

do not depend on t and have the form 

mi 

n(z) = 4>i[z){I itQ + (1 - 6 lA )I tA S t (z) + h, s d s - 2 p{z + v )), (3.3) 



s=2 



where 



m,i = i — 1, i = 1, . . . , I; rrii = I — 1, i = / + 1, . . . , N; (3-4) 

= t • (3-5) 

a\z — rj) l / z 

a(z + r\)a\z — (i — \)r\) 
In a neighbourhood of z = —r\ the function R(k, z) can be represented in the form 

l N 

R{k, z) = ]J(k + (z + r))- 1/2 hi{z + n)) ]J(k + (z + r,) 1/2 hi(z + rf)), (3.7) 

i=l 

where hi(z) are regular functions of z in a neighbourhood of z = 0. 

Proof. Due to ( 2.10| ) the matrix elements Lij are double-periodic functions of z. Therefore, we can 
consider them as single-valued functions on the original elliptic curve To with a branch cut between 
t he points z = ±77. First of all let us show that the coefficients Ti(z) of the characteristic poly nomial 
(3.1) are meromorphic functions on the Riemann surface To of the function E(z) defined by ( |2.9| ). (This 
means that are double-valued functions of z with square root branching at the points z = ±rj.) 



This fact can be immediately seen from the "gauge equivalent" form of L(t, z): 



L(t, z) = G(t, z)L(t, z)G-\t, z), Gij = 5 t 



a{z + rj) 



Xi(t)/2 V 



(3.8) 



where the matrix elements of L(t, z) have square root branching at the points z = ±77. Moreover, writing 
them explicitly, 

I. u z s _ ( b t a j) a(z + x l - gj) 

%J ' [<j{z - T])a(z + 77)] 1 / 2 a(x l - xj - 77) ' 

we conclude that rn(z) are single-valued meromorphic functions of z (i.e., elliptic functions). Further, 
T2i+\{z) are meromorphic functions on Tq, which arc odd with respect to the involution fo : To — » To, 
interchanging sheets of the covering To — ■> To (this involution corresponds to changing sign of the square 
root: E(z) —* —E(z)). Thus the curve T is invariant with respect to the involution 

f:fi — >t, f(k,E)\ — y(-k,-E), (3.10) 

which covers the involution fo- 

Let us consider the factor-curve: 

r:={f/f}. (3.11) 
This curve is a AT-fold ramified covering of the initial elliptic curve, 

r^r . (3.12) 

It can be defined by the equation 

N 

R(K,z) = K N + ^2r i (z)K n - 1 = 0, (3.13) 
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where 



K = k 



<t(z + 1]) 



1/2 



Ri(z) = n(z) 



a(z - rj) 



_o(z + T]) 



'72 



(3.14) 



Let us give some more comments on this statement. The coefficients Rj(z) of the equation ( 3.1 3j ) are 
meromorphic functions of the complex variable z obeying the following monodromy properties: 



Rj(z + 2lu u ) = R ] {z)e- 2 ' K ^^ . 



(3.15) 



Eq. (3.13) defines a Riemann surface f , which is a iV-fold covering of the complex plane. Due to ( 3.1 5| ) 
this surface is invariant with respect to the transformations 



2cj„, K 



(3.16) 



The corresponding factor-surface is an algebraic curve T, which is the covering of the elliptic curve with 
periods 2u a - 

Now we switch to ([3j]). This equality follows from the fact that the leading term of L(t, z) in a 
neighbourhood of the point z = —rj, 



Lij (t, z) 



Q>t°s) 



[<t(-2t,)(z + t,)]V'< 



+ 0((z + r 1 ) 1 l 2 ), 



has rank I. The corresponding N — Z-dimcnsional subspace of eigenvectors C = (ci, . 
eigenvalue is defined by the equations 



JY 



c i a 3' a ~ ^' 



a = 1, 



(3.17) 
, cat) with zero 

(3.18) 



Let us turn to determination of the coefficients in the characteristic equation (3.1). Since the matrix 
elements of L(t, z) have simple poles at the point z = n, then in case of general position the function 
ri{z) has pole of order i at this point. It follows from (3.7) that at the point z = —r\ the function r,;, 
i = 1, . . . , I, has pole of order i. If i = I + 1, . . . , 21 it has pole of order 21 — i. At last, if i = 21 + 1, . . . , N 
it has zero of order i — 21, i.e., 



n(z) = {z + ri) t/2 Pl (z + ri) 7 i 
n{z) = {z + riY/ 2 - l p l {z + j 1 ), i = 



l + l,...,N, 



(3.19) 
(3.20) 



where pi{z) are r egu lar functions. All these properties of the functio ns r j(z) allow one to represent them 
in the form ([D^, |3.4| ). (Note that the function (f>i(z) defined by eq. (|3.5| ) has pole of order i at the point 
z = r] and zero of order i — 2 at the point z = —rj.) ■ 



Imp ortant remark. It is necessary to emphasi ze t hat (3.7) implies that the characteristic equation 
(3.1) defines a singular algebraic curve. Indeed, (3/7) implies that (N — I) sheets of the corresponding 
ramified covering intersect at the point (z = —rj, k = 0). We keep the same notation T for the algebraic 
curve with the resolved singularity at this point. 



The coefficients Ii tS in ( |3.3| ) are integrals of motion. The total number of them is equal to Nl— 1(1— 1)/2 
which is exactly half the dimension of the reduced phase space. It follows from the results of the next 
section that they are independent. 



Lemma 3.1 In general position genus g of the spectral curve V (defined by eq. ( 3.1!% )) is equal to 
Nl- 1(1 + l)/2+l. 

Proof. First let us determine genus g of the curve T defined by eq. ( |3.l| ). By the Riemann-Hurwitz 
formula g we have 2g — 2 = 2N + is, where v is the number of branch points of T over Tq, i.e. , the number 
of values of z for which R(k, z) = has a double root. This is equal to the number of zeros of dkR(k, z) 
on the surface i?(fc, z) — outside the points located above the point z — —rj (due to the singularity 
of the initial curve mentioned above). The function dkR(k, z) has poles of order N — 1 above the point 
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z = rj It has also poles of the same order in I points located above the point z = — r\ that correspond to 
first I factors in ( |3.7| ). In the other N — / points above z — —rj it has zeros of order (N — l)(N — 21 — 1). 
Therefore, v = AIN — 21(1 + 1). The curve F is the two-fold branched covering of the spectral curve T, 
the number of branch points being equal to 27V (they are located above z = if]). The Riemann-Hurwitz 
formula gives then the relation 2g — 2 = 2(2g — 2) + 2N, which proves the lemma. ■ 



The characteristic equation (3.1) allows us to define two sets of distinguished points on the spectral 



curve (3.13). It follows from the factorization of R(k,z) (3.7) that the function k has poles at I points 



lying on the different sheets over the point z — —rj (they correspond to the first / factors in (3.7)). Let 



us denote them by P , i — Since a meromorphic function has as many zeros as poles, we 



conclude from ( 3.14 ) and (3.7) that the function k on the unreduced spectral curve T has 21 zeros which 
do not lie above the point z — —rj. These zeros correspond to I points P^~ , i = 1, . . . , I, on the spectral 
curve T: 

k(Pr) = 0. (3.21) 

In general position there is such a point Pf above each zero z~ of the function r^(z) different from its 
apparent zero z — —rj: 

(N-l)/2( , \ 2 

r N (z) = I„ i0 __J£^W J] o{z - z-). (3.22) 
In the abelian case (I = 1) the second marked point Pf~ lies above the point z = (N — 1)?/. 



Theorem 3.2 The components ^ a (x,t, P) of the solution ^(x,t,P) to the equation (2.1i) are defined 
on the N-fold covering T of the initial elliptic curve cutted between the points P^ and P^ , i = 1, . . . , I. 
Outside these branch cuts they are meromorphic. For general initial conditions the curve T is smooth. 
Its genus equals g = Nl — + 1 and ty a have (g — 1) poles 71, ... , J g —i which do not depend on the 

variables x,t. In a neighbourhood of P^ ', i = 1, . . . , I the function ^ a has the form 

00 

t, P) = [xt + E * a s 1 ^ + ilY)^ + r,)- 1 )''*^*^- 1 * *i(0, 0, P), (3.23) 

8 = 1 

where Xq 1 are constants independent ofx,t, and m are non-zero eigenvalues of the matrix (bfaj). In a 
neighbourhood of P~ the function ty a has the form 

00 

^ a (x 1 t,P)^(z-z-r^ (Y^xTM^-h f) *i(0,0,P) (3.24) 

s=0 



(z i are projections of the points P. to the initial elliptic curve; they are defined by eq. (3.2i)). The 



boundary values of the function Sfr a at the opposite sides of the cuts are connected by the relation 

*(+) = ^e 2mx/r >. (3.25) 

Proof. We start with analitic properties of the eigenvectors of the Lax matrix. 

By F* denote the curve T with cuts between the pre-images P^ of the point z = —rj and the 
pre-images Q~ of the point z = rj, i = l,...,N. For a generic point P of the curve F, i.e., for 



the pair (k,z) = P, which satisfies the equation (3.1), there exists a unique eigenvector C(0, P) of 



the matrix L(0, z) normalized by the condition Ci(0, P) = 1. All other components Ci(0,P) are given 
by Aj(0, P)/Ai(0, P), where Aj(0, P) are suitable minors of the matrix kl — L(0,z). Thus they are 
meromorphic on F*. The poles of Ci(0, P) are the zeros on T* of the first principal minor 

Ai(0,P) = det(fcJ y - Lij(0,z)) = 0, i,j > 1. (3.26) 

Therefore, these poles only depend on the initial data. 

Lemma 3.2 The coordinates Cj(0, P) of the eigenvector C(0, P) are meromorphic functions onT*. The 
boundary values of the functions Cj(0, P) at opposite sides of the cuts satisfy the relation 

c+ = c - e ™(* 3 (o)-*i(o))/^ (32 7 ) 
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In a neighbourhood of the point P^~ the junctions Cj(0, P) have the form 

Cj (Q,P) = (cf +) {0) + O{z + r 1 ))(z + T 1 )^^- x ^^, 
where cj* (i) are eigenvalues of the residue of the matrix L(0, z) at z = —rj, i.e., 

N 
3=1 

In a neighbourhood of the point Q~ the functions Cj(0,P) have the form 

Cj (Q, P) = (cf'-\0) + 0(z - r})){z - v )(^(o)-^(o))/2^ 

where Cj \t) are eigenvectors of the residue of the matrix L(0, z) at z = rj. 

The proof follows from the representation of (7(0, P) in the form (see (|3.8])): 

C(0,P) = G(0,z)C(0,P), 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



where G(t,z) is defined in (3.8), and (7(0, P) is the eigenvector of the matrix L(0, z). Matrix elements 
of L(0, z) ar e ana litic on the cuts betw een z = ± rj. Hence C(0,P) has no discontinuity on the cuts. 
This proves ( |3.27 ). The equations ( 3.28 ) and ( 3.3C ) are direct consequences of the fact that L(0, z) has 
simple poles at the points z — ±77. 

Remark. The vector (7(0, P) is invariant under the involution ( |3.1C| ) (that's why its argument is a point 
P of the spectral cur ve T rather than a point P € T). However, this notation is somewhat misleading 
since both factors in (3.31) are multi-valued on T, and only their product is well-defined. 

Lemma 3.3 The poles of C(0,P) are invariant under the involution f. The number of them is equal 
to 2NI -1(1+1). 

To prove the lemma we use the following standard argument. Consider the function of the complex 
variable z defined by: 

F( z ) — (Det |ci(0, Mj)|) 2 , 

where Mj, j — 1, . . . , N, are the points above z. It is well defined as a function of z since it does not 
depend on the ordering of the Mj's. The analitic properties of Cj allow us to represent F in the form 



F(z) = F(z) 



a[z + 77) 



(3.32) 



where F is a meromorphic function. This means that F has as many zeros as poles. The number of its 
poles is twice the number of zeros of the vector (7(0, P) whereas the number of zeros of F is equal to 
the number of branch points v of the covering T over Tq ( |3.l| ). In the proof of Lemma 3.1 we showed 
that v = 4NI — 21(1 + 1). The invariance of poles of (7(0, P) under the involution f follows from the 
f-invariance of eq. (3.26) which determines positions of the poles. This completes the proof. 

Let 71, . . . , 7g-i be the points of the spectral curve T whose pre-images are poles of (7(0, P). Note 
that if r is smooth, then g = Nl — Z(Z + l)/2 + 1 coincides with its genus. 



Let (7(t, P) be the vector obtained from (7(0, P) by the time evolution according to eq. (2.35). 

Lemma 3.4 The coordinates Cj(t,P) of the vector C(t,P) are meromorphic on T*. Their poles are 
located above the points 7i,...,7 g _i and do not depend on t. The boundary values of Cj(t,P) at 
opposite sides of the cuts satisfy the relation 



Ti(xj(t)-xi(0))/j7 



(3.33) 
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In a neighbourhood of P^~ the functions Cj(t,P) have the form 

Cj (Q, P) = (cf- +) {t) + 0(z + n))(z + r) ){MO)-*i<t))/*'i exp( Kl (z + V yH), 



(3.34) 



where Ki and c^' + \t) are defined in ( 3.2f\ ). In a neighbourhood of Q i the functions Cj(t,P) have the 
form 

c 3 (t,P) = [cf ~ ] {t) + 0(z - V ))(z - vY^^-^W 271 . (3.35) 



Proof. The fundamental matrix of solutions S(t, z) to the equation 

(d t + M(t, z))S{t, z) = 0, S{0,z) = l, 



(3.36) 



is a holomorphic function of the variable z outside the cut connecting the points z = ±77. We have 
L{t,z) = S(t,z)L(0,z)S- 1 (t,z). Therefore, the vector C(t, z) equals C(t,z) = S(t, z)C{0, z) hence it 
has the same poles as C(0, P). 



Let us consider the vector C(t, P) such that 

C(t,P) = G(t,z)C(t,P), 



(3.37) 



where G(t, z) is the same diagonal matrix as in (3.8). This vector is an eigenvector of the matrix L(t, z) 
and satisfies the equation 



(d t + M(t, z))C(t, P) = 0, M = G~ l d t G + G^MG. 



(3.38) 



The matrix elements of M are analitic at the cuts between z = ±77. Thus G(t, P) is analitic at the cuts 
on r. Therefore, the multi-valuedness of C(t, P) is fully caused by the multi-valuedness of G(t, z). This 
proves eq. (|3|). Eq. (HI) follows rom the analiticity of M at the point z = 7/. In a neighbourhood 
of z = —77 we have: 

M l3 {t,z) = ^- ) +0((z + r 1 f). (3.39) 
Therefore, in a neighbourhood of P^~ it holds 

d t c(t, P) = (m(t, z) + o(z a ))c(t, P), 

where 



IH(t,z) = ^(z + n) 1 +0(1), 



(3.40) 



are eigenvalues of the matrix M. This proves eq. (3.34) 



Let us now continue the proof. By the initial definition, 

N 

tf(z, t, P) = J2 s 3 (*, P)Hx- Xj (t), z)k x '\ Sj (t, P) - Cj (t, P)a 3 (t) , 



The solutions ^ of the linear problem ( |2.19| ) are defined on the curve T. In order to show that $ is well 
defined on the spectral curve T we use the equality 



Cj (t,P)$(x- Xj (t) , z)k x / r ' = 5j (t, P) 



a(z + x - xj + 77) 
a(z + r/)a(x ~ Xj) 



- rj) 
a(z + n) 



1/2I x h 



(3.41) 



Recall that components of the vector C are even with respect to the involution f ( 3 . 1 C| ) . The factor 



K(P) = k 



a(z - 77) 



-1 1/2 



o(z + rj) 



(3.42) 



is f-invariant too. Thus ^(x, t, P) is well defined on the spectral curve T. At the same time we see that 
poles of ^{x, t, P) coincide with poles of (7(0, P), i.e., they are located at the points 71, ... , J g -i- 
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Note that K(P) is a multi- valued meromorphic function on T with zeros and poles (which are 
nevertheless well defined) at the points P~ and , i — 1, respectively. Therefore, cutting T 



between P^, i — we can choose the branch of the third factor in ( 3.42| ) in such a way that '5 



becomes single-valued outside these cuts, and its boundary values at the sides of the cuts satisfy the 
relation (3.25). Consider now the behaviour of ^ in neighbourhoods of P^~ . In a neighbourhood of 
z = —r\ we have: 



a(z - 



■v) 



1 



a(z + T])a(x) 
Hence in a neighbourhood of P^~ it holds 



z + rj 



O(l) 




°jz - rj) 
a(z + Vj) 



1/2 



c/r) 



(3.43) 



(3.44) 



where fcj(z) is the branch of k(P) de fined by i-th factor of the equality (3.7). For i > I the product of 
the se cond and the third factors in ( 3.44 ) is regular in a neighbourhood of . Since the eigenvalues 
Ki in ( 3.29 ) equal zero for i > I, then the first factor is regular in a neighbourhood of P^,i > I, too. 



Therefore, the functions ^ a are regular at these points. Similar arguments for i — 1, . . . ,1, prove the 
equality ( 3.23 ). Note also that Vl/i(0, 0, P) in ( [3.23] ) has simple poles at the points i = 1, . . . , I. 



Consider now the term \q 1 m (3.23). By construction, it does not depend on x. Substituting the 
series ( Pf ) into ( Elf ), we see that it does not depend on t as well. ■ 



The following theorem can be proved by the same arguments. 



Theorem 3.3 The components ^> + ' a (x, t, P) of the solution 1 $> + {x,t, P) to the equation (2.2C) are de- 
fined on N-fold covering T of the initial elliptic curve with branch cuts between the points P^ and 
P^~ , i = 1, . . . , I. Outside these cuts they are meromorphic. In general position, the curve T is a smooth 
algebraic curve. Its genus equals g — Nl — l SktU _|_ \ an d /j fls (g _ X) poles 7^, . . . , T^Lu which do 

not depend on x, t. In a neighbourhood of P^~ , i = 1, . . . , I the function <f+> Q has the form 

00 

*+< a (x, t, P) = ( x p m + Xp ai (x, t)(z + ?7) s )Mz + t,)-!)-*'^-**^- 1 * V+'^O, 0, P), (3.45) 

s=l 

where Xq 1 are constants independent of x,t. In a neighbourhood of P^ the function ^+' Q has the form 



V + ' a (x,t,P) = (2 



y -\-x/r) 



(J2 * + ' x (0,0,P). 



(3.46) 



s=0 



The boundary values vp+' a; ( ± ) of the function at opposite sides of the cuts are connected by the 

relation 

$+.«;(+) = ^+>«:(-) e - 2 'r«/'7_ (3.47) 



Remark. Theorem 3.2 shows, in particular, that the solution ^> to the equation (2. IE) is (up to 
normalization) the Baker- Akhiezer function. In the next section we show that this function is uniquely 
defined by the curve T, the poles j s , the matrix xo an d the value xi(0). All these quantities are defined 
by initial Cauchy data and do not depend on t. However, it is necessary to emphasize that part of them 
depend on the choice of the normalization point toj that we have chosen as (to = 0). Any initial data 
{xi,ii,ai,bf\(bf ,ai) — ii} define the matrix L by means of eq. (2.32). The characteristic equation 
(|3.lD defines the curve T. Eq. (3.26) defines a set of g — 1 points 7 S on T. Therefore, there exists a map 



{xi,Xi,ai,bf\(bf,ai) = ±i} 1 — > {T, D 6 J(T)}, 
D = Y J A{ls)+xiU {0 \ 

s=l 



(3.48) 
(3.49) 



where A : T — ► J(T) is the Abel map and is a vector dep en ding on T only (see ( [4.13] )). The 
coefficients of the equation (3.1) are integrals of the system ( 1.21 )-( 1.23 ). As we shall see in the next 
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section the second part of the data ( 3.4S| ) define angle-type variables, i.e., the vector D(t) evolves linearly, 
D(t) = D(t ) + (t- t )C/ (+) , if a point of the phase space evolves according to eqs. ( |l.2l| )-( [l.23| ). These 
equations have the obvious symmetries: 



aiX ->. w- 1 ^, b+W, 



(3.50) 



where qt are constants and W is an arbitrary constant matrix. In the next section we prove that the 



data r, D uniquely define a point of the phase space up the symmetry transformations (3.5C) 



4 Finite-gap solutions of the non-abelian Toda chain 

Finite-gap solutions of the non-abelian Toda chain were found in the paper |n| by one of the authors. To 
construct the inverse spectral transformation for the spin generalizations of the Ruijsenaars-Schneider 
model, we recall the main points of this theory. Since we are working with a continuous variable x rather 
than the discrete variable n, some minor modifications of the construction are necessary. 

Theorem 4.1 Let T be a smooth algebraic curve of genus g with fixed local coordinates Wj t ±(P) in 
neighbourhoods of 21 points Pp, Wj } ±(Pj ) = 0, j = 1, . . . , I and with fixed cuts between the points Pj . 
Then for each set of g + I — 1 points 71, . . . , 7 9 +;-i in general position there exists a unique function 
ip a (%, T, P), a = 1 . . . , I, T = {tij-±, i = 1, . . . , 00; j = 1, . . . , 1} such that 

1°. The function ip a of the variable P 6 T is meromorphic outside the cuts and has at most simple 
poles at the points 7 S (if all of them are distinct); 

2°. The boundary values ipa of this function at opposite sides of the cuts satisfy the relation 

^+\x,T,P)=^-\x,T,P)e 2 ^. (4.1) 

3°. In a neighbourhood of the point Pp it has the form 

i> a (x, T, P) = wj x ± /v f£ £ s aj ' ;± (*, ZX±J exp (f^ /r,/ : /, /;: j , w jt± = w^ ± (P), (4.2) 

^ + (x,T)^S ai . (4.3) 

The proof of theorems of this kind, as well as the explicit formula for ij) a in terms of Riemann theta- 
functions, are standard in the finite-gap integration theory. We use the notation of the paper | p9| . 

It follows from the Riemann- Roch theorem that for each divisor D = 71 + • • • + 7 s +i-i in general 
position there exists a unique meromorphic function h a (P) such that the divisor of its poles coincides 
with D and such that 

h a (P+) = S aj . (4.4) 

The basis of cycles a°, b® on T with canonical intersection matrix being fixed, this function may be 
written as follows: 



where the Riemann theta-function 6{z\, . . . , z g ) = 6{zi, . . . , z g \B) is defined by the matrix B = (Bo.) of 
periods of holomorphic differentials on T; A(P) is the Abel map: J(F); 

9-1 

Rj = -JC-A(Pf)-J2A( 7s ), j = l,...,l (4.6) 

s=l 

9-1 

S i = -lC-A( lg _ 1+i )- ^A( 7s ), (4.7) 

s=l 
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and 

g+l-l I 

Z a = Z - A(P+), Z Q = -K- ]T A( 7i ) + A(P+), (4.8) 

i=l 3=1 

where K, is the vector of Riemann's constants (the proofs can be found in |fL9| ), 

Let rff2^ lJ;± ^ be the unique meromorphic differential holomorphic on T outside the point P^, j — 
1, . . . , I, which has the form 

dfiW;*) = d(wjX + 0(wj.,±)), (4.9) 
near the point Pf 1 and normalized by the conditions 



It defines a vector [/( l J ;± ) with coordinates 



MW,±) =J-I dil^l (4.11) 
2tti J b o 

Further, define a differential dn^°\ which is holomorphic outside the points Pp, has the form 

dfi(o) = ±*£i± + 0(1)^. ±) (4.12) 

near these points, and has zero a-periods. It defines a vector with coordinates 

It follows from Riemann's bilinear relations that 

i 

C/(°) = ^- 1 ^(A(Pr)-^(P+)). (4.14) 

3=1 

Theorem 4.2 TTie components i]j a of the Baker- Akhiezer function ij)(x, T, P) are given by the formula 
4 -u (r) e{A{P) + Ui ° )x ± U (A) tA + Z a )9(Z ) ( , n (°) (P)+E a urf^P,) , . 

flW(P)= /" dfi^ A=(i,i;±). (4.16) 

Remark. Note that since the abelian integral f^°) has logarithmic singularities at the marked points, 
one may define a single- valued branch of ip only after cutting the curve T between the points P- . 



Let us now define the dual Baker- Akhiezer function. For each set of g + l — 1 points in general position 

)i 
3 



there exists a unique differential dft holomorphic outside the points Pf such that it has simple poles at 
these points with residues ±1, i.e., 



dCl = ±-^± + 0(l)dwj ± (4.17) 
and such that it equals zero at the points "f s : 

dQ(j s )=0. (4.18) 
Becides 7 S this differential has g + l — 1 other zeros which we denote by 7+. 
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The dual Baker- Akhiezer function is the unique function ip + (x, T, P) with components tp + ' a (x, t, P) 
such that 

1°. The function ip + < a of the variable P G T is meromorphic outside the cuts and has at most simple 
poles at the points 7+ (if all of them are distinct); 

2°. The boundary values ■0+ ;Q '( ± ) of this function at opposite sides of the cuts satisfy the relation 

r /,+;a,(.+) fa T) p) = ^+;«,(-) j, p) e -2«x/r,. ( 4 19 ) 



3°. In a neighbourhood of f^. it has the form 

^+< Q (*,T,P) = W fX /V ■' J: -'■■!■■< ^ ( X "I..':'',-) > ( 4 - 20 ) 

^+( a: ,T) = 5^. (4.21) 

Let h+{P) be the function that has poles at the points of the dual divisor 7+, . . . , 7 fl l + ;_ 1 and nor- 
malized by /i+ (Pj~) = $aj- It can be written in the form (4.5) in which 7 S are replaced by 7+. It follows 
from the definition of dual divisors that 

9+2-1 9+2-1 1 

Y, A( ls )+ £ A( 1 t)=K + J2(A(P+) + A(Pp), (4.22) 

S—l S—l j — 1 

where Kq is the canonical class (i.e. the equivalence class of the divisor of zeros of a holomorphic 
differential). Thus the vector Zq in the formulas for /i + is connected to Z by the relation 

l 

Z + Z+ = -2K - K + Y(A(Pf) - A(Pr)) = -2K -K - Ef<%. (4.23) 

i=i 



Theorem 4.3 The components ip + (x,T, P) of the dual Baker- Akhiezer function are given by 
M h+(P] WP) U (0) * ~ Ea U {A) tA + Zj)6{Z+) {xam(P)+T 

v a( ) 9(A(p) + z+)e(u^x + Y,AU^t A ~z+) e ' l4 -' 4J 

where 

= -Z - 2K - K - Z+ = Z+-A(P+). (4.25) 



The above results are valid for any algebraic curve with two sets of marked points. Consider now 
the class of curves corresponding to the spin generalizations of the Ruijsenaars-Schneider model. 

Theorem 4.4 Let T be a smooth algebraic curve defined by an equation of the form 

N 

R(K,z) = K N + J2R l (z)K N - 1 = 0, (4.26) 

i=l 

where Rj(z) are meromorphic functions of z such that 

Rjiz + 2tu a ) = R j (z)e- 2iC ^ , i. (4.27) 

and holomorphic in the fundamental domain of the lattice with periods 2uj a outside the point z = —r). 
Let us assume that in a neighbourhood of z = —T) the polynomial R has the following factorization: 

l N 

R(K, z) = Y[(K + (z + r,)- x Hi(z + 77)) {K + (z + ry)^(z + 77)), (4.28) 

i=l i=2+l 
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where Hi{z) have no singularity at the point z = —X). 

Then the Baker- Akhiezer function ip corresponding to: (i) the curve T, which is the factor ofT with 
respect to the transformation group 

zi — >z + 2uj a , if i — ► Ke- 2 ^^, (4.29) 

(ii) local coordinates Wj.+ — (z + T])Hj(Q) near the poles Pj~,j = 1, • ■ ■ , I of the multivalued function 
K = K(P) and arbitrary local coordinates Wj t — near the zeros PJ of this function obeys the relation: 

iP(x + 2u a ,T,P)=cp a (P)ib(x,T,P). (4.30) 

where 

( p a (P)=K(P) 2u '"/ ri e c ^ a >. (4.31) 



The proof is easy. It follows from the monodromy properties ( 4.29 ) that the values of <p a (P) do not 
change under shifts of z by periods of the elliptic curve, i.e., it yields a well-defined function on T. Eq. 
(4.30) follows from the fact that its left and right hand sides have the same analitical properties. ■ 



Corollary 4.1 The Baker- Akhiezer function ip(x,T,P) with components tj) a (x,T, P) 
corresponding to the data of Theorem 4-4 can be written in the form 



o{z + ri) 



i=i 

Under these very conditions the dual Baker- Akhiezer function has the form 

m 

il>+(x, T,P)=J2 st(T, P)*(-x + Xi(T) - rj), z)k~'^. 



1/2 



(4.32) 



(4.33) 



Proof of the lemma is identical to the proof of the corresponding lemma from the paper J19| and it 
was already presented at the beginning of Sect. 2, because Theorem 4.4 implies that ip and tp + arc 



double-Bloch functions. (Note that from (4.23) and (4.24) it follows that ip + has poles at the points 
Xi - rj). ■ 

So far ti t a;± were arbitrary parameters, entering tp through the form of the essential singularity at 
Pt- Fix now some values of these parameters for i > 1, i.e., ti jQ; ± = t® a . ± , whereas for i = 1 we put 



h,a;± = t± +t 1 j 



(4.34) 



The Baker- Akhiezer function now depends on the variables (x, t-\-,t-). For the sake of brevity, we denote 
it by tp(x, t-|_, t_, P), skipping the dependence on the constants T°. 



Theorem 4.5 For any choice of the constants T° the Baker- Akhiezer function 
ij)(x,t+,t—, P) satisfies the equations 

d + tjj(x, t + , t-,P) = tp(x + rj, t+,t-,P) + v(x,t+,t-)ip(x,t+,t-,P), 

d-ib(x, t + ,t-,P) = c(x, t + , t-)ip(x — 7], t + , t-,P), d± — d/dt±, 

where 

v(x, t+,t-) = d+g(x, t + ,t^)g^ 1 (x, t+, i_) = ££(x, t + ,t-)-£(x + r), t+,t-), 
c(x,t+,t-) = g(x,t + ,t^)g~ 1 (x - 7?,* + ,i_) = d-^(x,t + ,t-), 
and the matrices g and £^ are defined by the coefficients of the expansion ({A ): 



(4.35) 
(4.36) 

(4.37) 
(4.38) 



g°"i(x,t + ,t-)=& j >-(x, t + ,i_), ti(x,U,t-) = {tf r ' + } 



(4.39) 
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Proof. Eq. (4.35) follows from the fact that the function 

d + 1p a (x,t+,t-, P) -1p a (x + T],t+,t-,P) 



has the same analitic properties as ip except the normalization 4.3). Thus we can write it as a linear 
combination of the components tftp with coefficients v a @ . The explicit form of v a @ follows from comparing 
the coefficients in the left and right hand sides of ( 4.35 ) at the points PJ . Hence we get the first equality 



in (4.37). On the other hand, we may expand ip near the points P- . Then we get the second equality 
in (4.37). Eqs. ( 4.36| ) and (4.3S) are proved in a similar way. ■ 



Corollary 4.2 The matrix function g n (t+,t_) = g(jirj + Xq, t+, t_), corresponding (according to the 
definitions of the Baker- Akhiezer functions ) to the Riemann surface T with fixed local coordinates near 
the marked points Pp and to the set of points j%, . . . ,g+l — 1 is a solution of the 2D Toda chain equations 



Remark. The dependence of g n on the variables tij, 
flows of the 2D Toda chain hierarchy. 



,oo, j = 1, 



, I, corresponds to higher 



Theorem 4.6 The dual Baker- Akhiezer function satisfies the equations 

- d+i> + (x,t + ,t-,P) = %l)+{x -r),t+,t-,P) + ip + (x,t+,t-,P)v(x,t+,t-), 
- cL V;+ (x, t+ , t- , P) = ip+ (x + V , t+ , t_ , P)c{x + n, t+ , t _ ) , 



(4.40) 
(4.41) 



where c(x,t+,t—), v(x,t + ,t-) are the same as in (4-35, 4-.3L). 



Proof. The same arguments as in the proof of Theorem 4.5 show that if> + obeys eqs. ( 4.4C| ), (4.41) with 
coefficients v + and c + given by 

c+(x,t+,t-) = [^-(x + V ,t + ,t^)}- 1 ^-(x,t + ,t^), (4.42) 

v+{x,t+,t-) = -lH+ ; -(x 1 t + ,t_)]- 1 d + t;+'-(x,t + ,t^, (4.43) 

where the matrix elements = {CcT'""'' - } are determined from ( 4.20| ). The coincidence of the coeffi- 
cients for ip and ip + is easily seen from the relation 



(4.44) 



which follows from the definition of the dual Baker- Akhiezer function. To prove (4.44), consider the 
differential ip a ip + Pd£l, where dfi is the same as in the definition of the dual divisor 7+. It is a meromorphic 



differential on T with the only poles at P- . Its residue at P- is 



res 



i+M+Pdn = 6 

a, 



Further, the residue of this differential at Pj is 



j 

Since the sum of the residues should be zero, then ( [4.44 ) follows from (4.45) and (4.46) 



(4.45) 



(4.46) 



(4-Sl, i-3t ) are elliptic functions in x. In general position they have the form 

N 



Theorem 4.7 Let the curve Y , the marked points Pp and local coordinates in their neighbourhoods be 
the same as in Theorem 4-4> then the corresponding algebraic- geometrical "potentials" v and c in eqs. 

(4.47) 

(4.48) 



v(x, T)=J2 a*(?>+(PMz - Xi(T)), 

i=l 

2(x,T) = d- (^(T) + J2'H{T)bf{T)C{x-x i (T) y j , 



where at, bj are some vectors, and Sq is a matrix functions that does not depend on x. 
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Proof. The potentials in eqs. ( 4.3 5| ) and (4.36) arc elliptic functions due to (4.30). It follows from ( 4.1 5| ) 



that the poles x = Xi (T) of the Baker- Akhiezer function correspond to solutions of the equation 

9(U (0) x + J2u iA h A + Z Q ) =0. (4.49) 



I t is clear from ( |4.8| ) that for corresponding solutions (xi(T), T) the first factor in the numerator of 
( 4.15 ) is zero at P = P+. At the points P^ , (3 ^ a the function h a {P) vanishes. Therefore, the 



residue ip^ t (T, P) of the function ipa(x, T, P) at x = Xi(T) (as a function of P) has the following analitic 
properties: 

1°. It is a meromorphic function on T outside the cuts between the points Pp and has the same poles 
as ip; 

2°. Its boundary values ip° a \ (T, P) at opposite sides of the cuts satisfy the relation 

^ 3 + \T,P) = <^ ) (T,P)e 2 ^( T )/"; (4.50) 
3°. In a neighbourhood of P. we have 

oc 

<,(T, P) = wj; ± ^ exp(£ wJp. d . ± )F^ a (w Jt± ), (4.51) 

s=l 

where a are regular in these neighbourhoods, and 

F+J0) = 0. (4.52) 

Hence we see that ip^ i has the same analitic properties as the Baker- Akhiezer function except the 
following. The regular factor in the expansions of this function near all the points Pj~ has the vanishing 
leading term. For general x,tA there is no such function. For the special values (x = Xi(T),T) such 
function ipio (T, P) does exist and is unique up to a constant (with respect to P) factor (it is uniquely 



defined in general case when Xi{T) is a simple root of the equation (4.49)). So we can represent ip a in 
the form 

Mx , T , P) = 0a(T) + 0((x - Xi (T))°). (4.53) 

x — Xi(T) 

It follows from the last equality that the residues Pi(T) of the matrix ^(x,T) with matrix elements 
^ ;+ (x,T) 

Zt(x,T)= P ^ T) 1 + ((x-x i (t))°) (4.54) 

X Xi\l ) 

have rank 1. This means that there exist vectors a^(T) and covectors bf(T) such that pi = ai(T)bf (T). 



We see from ( 4.3C ) that the matrix obeys the following monodromy properties: 

+ 2w,) - tffoT) + 2C(w,)r, (4.55) 



where r is a constant. It follows from these equations and (4.54) that £^ can be written in the form 

N 

= S (T) + ai(T)bi(T)((x - Xi (T)). (4.56) 
Combining this with the second equalities in ( [4.37j ) and (4.38), we get the desired result. ■ 



Remark. In the abelian case (I = 1) there is an equivalent expression for c(x) as a product of pole 
factors: 

c^T)=n g(x "^ (T) 2 r ??M ^r (r) ' ?7) 

fj[ o- 2 {x-Xi{T)) 
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(see (1.17)). Comparing the pole terms in (4.48) and (4.57), we get the relations 

a(xi - x k + rj)a(x l - x k - ?/) 



d+d-X i (T) = -a 2 (r ] )l[ 

k^i 



<J 2 (Xi - X k ) 



d+d^Xi(T) = -d + Xi(T)d-Xi(T)J2(V(xi ~ x k ) ~ V{x k - x { )). 

k^i 



(4.58) 
(4.59) 



It is easy to check that if the dynamics with respect to t± is given by the hamiltonians a(±rj)H± ( 1.20 ), 
then these relations follow from equations of motion. 

The connection between algebraic-geometrical potentials in eqs. ( 4.35Q , (4.36) corresponding to 
equivalent divisors is well known in the theory of finite-gap integration. Let D = *yi + • • • + 7 g +;-i 
and = j^' + • ■ • + 7^ ; _ 1 be two equivalent divisors. This means that there exist a meromorphic 
function h(P) on T such that D is the divisor of its poles and D^- 1 ' is the divisor of its zeros. 



Corollary 4.3 The algebraic-geometrical potentials v(x,T), c(x,T) and v^'(x,T), 

c^(x,T) corresponding to T, Pj ,wj t ±(P) and equivalent divisors D and D^ 1 ' are gauge equivalent: 



{x,T) = Hv(x,T)H- 



(x,T) = Hc(x,T)H-\ H aj = h{P+)5 aj 



(4.60) 



Corollary 4.4 In general position the curve T satisfies the conditions of Theorem 4-4 if and only if it is 
the spectral curve ( 3.1^ ) of the Lax matrix L defined by eq. (2. Si ), in which arbitrary constants 

and the vectors aj, b^ satisfy the conditions (1.25). 

It follows from Theorem 3.2 that the Baker-Akhiezer function ^5> a (x,t, P)/ 1 i>i(0,0, P) is related to 
the normalized Baker-Akhiezer function ip a (x,t, P) by the formula 



V a (x,t,P) a/3 , , „, 



(4.61) 



where ip(x,t,P) is the Baker-Akhiezer function defined in the beginning of this section. It corresponds 
to the following values of the parameters T = {£j,j;±}: 



tl,j;+=t, ti,j,±=0, (i,3,±) 0-,3,+), 



(4.62) 



The equality (4.61) yields the corollary: 

Corollary 4.5 Let ai(t) , bi(t) , Xi(t) be solutions to the equations of motion (1.21)-( 1.23\ ), then 

N 

J2<t)b+(t)V(x-x l (t))= X ov(x,t)xo\ (4.63) 



where v(x,t) = v(x, t+ = t, t_ = 0) is the potential corresponding (according to Theorem 4.5) to the 
normalized Baker-Akhiezer function %j)(x,t,P) constructed from the data (a curve with marked points) 
obeying the conditions of Theorem 4-4- 



Corollary 4.6 The map 

ai (t), bf{t),x 4 (t) .— > {r, [D]}, (4.64) 
where [D] is the equivalence class of the divisor D is an isomorphism up to the transformations ( 3.5(\ ). 

The curve T does not depend on time, while [D] depends on the choice of the initial point to = 0. The 
follwoing theorem shows that this dependence [D(io)] is linearized on the Jacobian. 
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Theorem 4.8 Let V be a curve that is defined by the equation ({.2h. ) and D = 71, . 
lg+1-1 be a set of points in general position. Then the formulas 



(4.65) 



= Qi 1 {t)h a (q ) 



i(uWxj(t) + u<-+k + z a ) 

0(Z a ) 



tl(Z a ) 



where 



1 1 

Q 2 S) = ~Y, h t(<*>) h M 



0(U^Xj(t) + U {+) t - Z a )8{U^ Xl (t) + U {+) t - Z+) 
9(Z a )8(Z+) 



(4.66) 
(4.67) 

(4.68) 



(go is an arbitrary point of the curve T) define solutions of the system fiTjTl ), ( 2. 25), (2.24)- Any solution 
of this system in general position may be obtained from ( 4.6t\ )-(4-6'i ) by means of the symmetries ( 3.5C ). 



Remark. If in eqs. ( 4.65| )-(4.67) the vector Zq is substituted by 

Z .— > Z + J2u (A) t A , 



(4.69) 



then the corresponding quantities Xj(T ), aj ( T), h (T), T = {t^}, depende on tA according to the 
higher commuting flows of the system ( |l.2l| )-(l .23 ). We would like to emphasize that the points 
enter symmetrically. This means that the dependence of Xi(T), ai(T), bi(T) on the variable t- = 
I 1 Si=i hj;- ^ s described by the same equations as that for t = t+. 



5 Difference analogs of Lame operators 

Consider the operator So given by eq. ( p. .9] ) for integer I. First of all we note that due to the obvious 
symmetry —I <-> £ — 1 it is enough to consider only positive values of i. In what follows we imply that 
I G Z+. The finite-gap property of the operator So means that the Bloch solutions to the equation 

(S Q f)(x)=ef(x) (5.1) 

are parametrized by points of a hyperelliptic curve of genus 11. 
Any solution f{x) of eq. flj.ip may be represented in the form 



fix) = *{x) (ei(v/2)r x/v II hi? ~ jri), (5.2) 

3=1 

where ^(x) satisfies the equation 

(5b*) (a:) = *(« + V) + ce(x)V(x -rj) = s^(x), (5.3) 

c,(x)= 9 ?(,/2)^±^(f^±M. (5.4) 

This transformation sends Bloch solutions of the first equation to Bloch solutions of the second one. To 
construct them explicitly, we use the ansatz similar to the one for the linear equation ( |2.2| ): 

l 

* = J2 S J ( z ' fc > £ )^( x - 3V, *W h , (5.5) 



24 



where 



$t(x,z) = 



\(z + x + N-q) 
{ {z + Nr))8x(x) 



x/2r) 



N = 



1) 



(note that $1 coincides with $(x, z) given by (2.6) and c\[x) coincides with c(x — r/) from 
constant factor). 

The function <$>t(x, z) is double-periodic in z: 

$e(x,z + 2u a ) = $t{x,z). 



(5.6) 
up to a 

(5.7) 



(In this section we take the periods to be 2lo\ = 1, 2ll>2 = t.) For values of x such that x/2r) is a 
half-integer nu mber the function $^ is single-valued on the Riemann surface Tq of the function E(z) 
defined by (2.9). For general values of x one can define a single-valued branch of <&t(x, z) by cutting the 
elliptic curve To between the points z = ±77. 

As a function of x, <&e(x, z) is a double-Bloch function: 

$e(x + 2uo a , z) = T(')(z)$ f (x, z), (5.8) 

where 



i( z - v) 
t(z + v) 



T 2 (£) (z) = exp(-27ri(z + Nr/)) 



i/2n 



y{z - v) 



(5.9) 



(5.10) 



h(z + v), 

In the fundamental domain of the lattice defined by 2uj a the function <$>i(x, z) has a unique pole at the 
point x = 0: 



1 



?[(0)x 



0(1). 



(5.11) 



Substituting (5.5) into (5.3) and computing the residues at the points z — jrj, j = 0, .. .,£, we get 
I + 1 linear equations 



3=1 



0, i = 0,. 



for £ unknown parameters Sj = Sj(z,k,e). Matrix elements Li j of this system are: 

L Q . l = k + h^,X-2f 1 ,z)k- 1 , L ,j ■ = h$t(-(j + l)ri,z)k- 1 , j = 2,...J; 

Li,i =-e- h<f> e (-r), z)^ 1 , L h2 =k — h$ t {-2v, L ltj = -h$t(-jr}, z)k~ x , j > 2; 

= 8ij+iCik~ - eSi t j + 6ij-ik, i > 1, 



where 



h = 0i(O)res :c= oc,( a ;) = 6 jM^6 l {£ri)9 1 {{l + 



(5.12) 

(5.13) 

(5.14) 
(5.15) 

(5.16) 
(5.17) 



The overdetermined system (5.12) has non-trivial solutions if and only if rank of the rectangular matrix 
Li_j is less than £. By and L^ 1 ' denote I x £ matrices obtained from L by deleting the rows with 
i = and i = 1, respectively. Then the set of parameters z, k,e for which eq. (5.3) has solutions of the 
form ( 5.12j ) is defined by the system of two equations: 



deti w = R {V, { Zl k,e) = 0, £ = 0,1. (5.18) 
Expanding the determinants with respect to the upper row, we get that R^ have the form: 



R (0 \z, k, e) = rf\s) + £ k^^Hv, ^(e), 



(5.19) 
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(z, k, e) = kh-i(e) + k-^ti-U + zyfl^e), 



(5.20) 



1 and £ — j respectively. 



where f^_i, r^-, r^. are polynomials in £ of degrees 

These equations define an algebraic curve T realized as £{£+ l)/2-fold ramified covering of the genus 2 
curve To on which the functions $£ (— jr/, 2) are single- valued. This curve enjoyes the obvious symmetry 

(z, fc, s) 1 — ► (z,-fc,-e). (5.21) 

This is a direct corollary of the following general property of eq. fl5.3| ). Let ^(x) be a solution of ( |5.3| ) 
with an eigenvelue e, then exp(7ria;/ry) is a solution of the same equation with the eigenvalue — e. 
This transformation corresponds to the change of sign for k. 



Note that the transformation k — > — fc in (5.5) and the simultaneous interchanging of sheets of E(z) 
leaves the function ^(x,z,k) invariant. Therefore, T may be considered as a £{£ + l)-fold ramified 
covering of IV 

Let us show that this curve is invariant with respect to another involution, 



(z, k, e) (-z, fc- 1 ^( ?7 /2), e), 



as well. 



Let ^(x) be a solution of eq. (5.3). Then the function 

*(af) = *(-a;)A(x), 



where 



A(x)=e 1 ( ?7 /2) 2 ^n 



2s/„ TT + 



=1 01 ~ ^) 



(5.22) 



(5.23) 



(5.24) 



is also a solution of the same equation. Moreover, if ^ is a double-Bloch solution then iff is also a 
double-Bloch solution. It is easy to check that if Bloch multipliers of ^ are parametrized by the pair 
z, k, then the Bloch multipliers of correspond to the pair —z, k~ 1 9\{q/2). 

The variable e being considered as a function e{P) on the curve T (P G T) is a meromorphic function. 
It can not take any value more than twice because for a given value of e the second order difference 



equation (5.3) has at most two different Bloch solutions. The involution (5.22) is not trivial, so the 
function s does take generic value two times. Therefore, the algebraic curve T is a hyperelliptic curve of 



finite genus g. Due to the symmetry (5.21) it can be represented in the form 



9+1 



(5.25) 



The involution (5.22) is the hyperelliptic involution corresponding to the interchanging of the sheets of 
the ramified covering J5.25 ). 



Now we are going to prove that g = 2£. From (5.25) it follows that there are 2g + 2 fixed points of the 
hyperelliptic involution. On the other hand, the number of fixed points of the hyperelliptic involution 
(5.22) is equal to the number of preimages of the second-order points u a £ To, a = 0, ... ,3, (that are 
fixed points of the involution z — > —z on Tq) on T such that the corresponding value of k is equal to 
±#1(77/2). Explicitly, from now on we adopt the following notation 



u = 0, wi = 1/2, wa = (1 + r)/2, lo 3 = r/2 . 



(5.26) 



For values of e corresponding to the fixed points there is only one Bloch solution. This implies that 
the corresponding solution has a definite parity with respect to the involution ( 5.23| ), i.e. 



= uty(-x)A(x), v = ±l. 



(5.27) 



1 this notation differs from that adopted in Sect. 2. 
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We are going to prove that v = (— l) e 



The equality (5.27) for x — rj implies that 

Sl = ^k- 1 9(-T } )(-iy- l ^ll^-9 1 (£r ] )d 1 ((£ + l)n), 



(5.28) 



Comparing this equality with ( p.3[ ) (taken at x — 0), we see that if si ^ 0, then v = (— 1) . Otherwise, 
si = and ^(-rj) = 0. 

From ( |5.3| ) it follows that coefficients si and S2 in the representation ( |5.5| ) of a double-Bloch solution 
can not equal zero simultaneously. Indeed, let j be a minimal index such that Sj ^ 0. Suppose j > 2, 
then the left hand side of (5.3) has a pole at the point z = (j — l)r/ but the right hand side has no pole 
at this point. 



T herefo re, if s% = 0, then S2 ^ 0. In that case eq. (5.2) implies that ^(O) ^ 0. At the same time 
from ( 5.27 ) (taken at x = 0) it follows that 

*(0) = (-iyW(O). (5.29) 

Therefore, v = (— 1) . 

For z = u a , k = ( — l) ct+1 6 l i (77/2) the representation (|5.5|) of the double-Bloch functions is equivalent 



to 



where 



= (6 1 (rj/2))^ exp(7rz(5 2 , + 5 3<a )x) f[ + % 



E 



•'• i = Wo 



(5.30) 



(5.31) 



Lemma 5.1 TTie hyperelliptic involution ofT has 2d fixed points, where d equals sum of dimensions of 
the functional subspaces consisting of functions that have the form (5. SL) and satisfy the relation (5.21) 
with v — (— l) e . 



Proof. As it was shown above, for each p air o f fixed points of the hyperelliptic involution of T that are 
invariant with respect to th e in volution ( 5.21 ) (co rresp onding to the ch ange of sign for k) there exists 
the unique solution to eq. (^3) having the form ( 5.30| ) and satisfying ( 5.27 ) with v = (— l) 1 . On the 
other hand, the space of such functions is invariant with respect to the operator Sq. Indeed, the equality 
( 5.28 ) with v = (— l) e (that is a corrolary of ( 5.27 )) implies that So^> has no pole at z = 0. At t he same 
time iSo commutes with the linear operator (5.23). Therefore, the number of solutions of (5.3) having 
the form ( 5.3C| ) and satisfying ( [5.27 ) is equal to the number of eigenvalues £.;, i = 1, . . . , d, of So on these 
finite-dimensional spaces, i.e., is equal to sum of their dimensions. ■ 



It is easy to see that dimension of the space defined in the last lemma for loq — is equal to (£ — l)/2 
if I is odd and 1/2 +1 if I is even. For the other three points of second order corresponding dimensions 
are equal to [l + l)/2 if i is odd and £/2 if £ is even. Therefore, the number of fixed points is A£ + 2 and 
we do prove that g = 2£. 



Lemma 5.2 The direct sum of spaces of functions having the form ( 5.31 ) and satisfying ( 5.2\ ) with 
v = (— l) e is invariant with respect to the operators S a , 



(S a y)(x)=H a 



9 a+1 (x - £n)_ nx + _ eKrj/2) 6 X (x-(£ + 1)t,)6 a+1 (-» - A?) _ ^ 



i(x-£n) 



i(x)6i(x - 77) 



(l) ^A+i(W2) 



(5.32) 
(5.33) 



which are gauge equivalent to the operators 

2 cf. the "sum rule" in |35^ 
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The proof is straightforward. Again, eq. (5.28) implies that S a ^> has no pole at x — 0. At the same time 
these operators commute with the transformaton (5.23) and keep invariant the set of Bloch multipliers 
corresponding to the spaces of functions having the form ( |5.30 ). Note that coefficients of S a , a ^ are 
not elliptic; that's why they do not preserve each space but only their direct sum. 

Remark. The invariant space for the operators S a described above coincides (after the gauge transfor- 
mation (5.2)) with the finite-dimensional representation space of the Sklyanin algebra found in More 
information on invariant functional subspaces for Sklyanin's operators is contained in the next section. 



Equations (5.18) represent k and e as multi- valued functions of z. Let us consider now analytical 
properties of the double-Bloch solutions ^(x,Q), Q = (z, k, e) £ T on the hyperelliptic curve T. 



Theorem 5.1 The Bloch solution ip(x,P) = ^(x,P)^> 1 (0, P) of eq. \5. 3\ ) is a meromorphic function 
on the genus 2£ hyperelliptic curve T defined by the equation 



21 + 1 



(5.34) 



(where E j are eigenvalues of Sq on the finite- dimensional invariant space of functions having the form 
( 5.3(\ ), ( 5.3\ )) outside a cut between the points P ± at infinity (that are preimages of e = oo on f). 
Outside the cut ip(x,P) has 2£ poles independent of x. They are invariant with respect to the involution 
of f 1 that covers the involution e — > — e (5.21). The boundary values of ip^'(x,P) at opposite sides of 
the cut are connected by the relation 

tf,W = ^-) e 2 ^ /r '. (5.35) 
In a neighbourhood o/P ± the function ip(x,P) has the form: 



\s=0 
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(5.36) 



Proof. Coefficients Sj in (5.5) are solutions of the linear system (5.12). Let us normalize them by the 
condition s± = 1. Then all other Sj are meromorphic functions on T. Therefore, we may conclude that 
^ (as a function of Q £ L) is well defined on T with cuts between zeros and poles of 



K = 



v(z - rj) 



1/2 



h{z + rj). 

At the edges of these cuts ^ has a singularity of the form 



k{z). (5.37) 



(5.38) 

(up to a factor of order 0(1)). From (|5.3| ) it directly follows that a function ^ with such type of 
singularity might be a solution of this equation if e = oo at the s ingul ar points. On the other hand, in 
neighbourhoods of these points we have K ~ e^ 1 . That proves ( |5.36|) . The proof that the number of 
poles of ip is equal to 2£ is standard. First of all note that poles of t/j do not depend on x because poles 
of Sj do not depend on x. Then we may consider the following meromorphic function of e: 



F(e)^^( Vl P 1 (e))-i;(r,,P 2 (e))f 



(5.39) 



where Pi(e) are two preimages of e on T. This function does not depend on the ordering of these points, 
i.e., F is a meromorphic function of e. It has double poles at projections of poles of tp, has a pole of 
order 2 at infinity and has simple zeros at the branch points. The numbers of zeros and poles of a 
meromorphic function are equal to each other. That implies that ip has g poles. ■ 



We would like to mention that the theorem could be proved by a direct consideration of analytic 
properties of 'J on T represented in the form given by (5.18). This alternative proof is very similar to the 
proof of Theorem 3.2. We skip it, but present some arguments explaining why the multi- valued function 
K has only one pole and one zero. 
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For I > 1 the function <J>^ (— jrj, z) has a pole and a zero of order j at the points z — ±77, respectively. 
Therefore, k has zeros and poles at all preimages of z = —r\ (resp., z — rj) on T. Hence K is regular at 
all these points. The function $^ has a simple pole at the point z = —Nrj. It turns out that at one of the 
preimages of this point the function s (as well as k) has a pole. The corresponding point is one of the 
infinities on T (representated it in the form (5.18)). The other point at infinity is one of the preimages 
of the point z = Nrj. 

Remark. Eqs. ( |1.9[ ) define one particular series of representations of the Sklyanin algebra found in || 
(called series a) there). In this series the operator So has the finite-gap property. It should be noted 
that in each of the other series there is again an operator having the finite-gap property (see (6.31) 
below). The other three operators in each series in general do not have this property. Trigonometric 
degenerations of these operators were considered in |32j ] . In this case So corresponds to solitonic solutions 
of the Toda chain. 



6 Representations of the Sklyanin algebra 

In this section we show how to construct representations of the Sklyanin algebra from "vacuum vectors" 
of the L-operator ( |l.4| ). The notions of vacuum vectors and the vacuum curve of an L-operator were 
introduced by one of the authors |?3| in analysis of the Yang-Baxter equation by methods of algebraic 
geometry. Let us recall the main definitions. 

To begin with, consider an arbitrary L-operator C with two-dimensional auxiliary space C 2 , i.e., an 
arbitrary In x In matrix represented as 2 x 2 matrix whose matrix elements are n x n matrices A, B, 
C, V: 

Z=( A r SV (6-1) 



C V f 

The operators A, B, C, T> act in a linear space TL = C which is called the quantum space of the L- 
operator. We emphasize that so far no conditions on C are implied. In particular, we do not impose the 
relation ( |l.5| ) and do not imply any specific parametrization of the matrix elements. 

Let us consider a vector X ®U £TL®C 2 (X ETL, U £ C 2 ) such that 

C(X ® U) = Y <g> V, (6.2) 

where Y £ TL, V £ C 2 are some vectors. The relation ( |6.2[ ) means that the undecomposable tensor 
X (g) U is transformed by C into another undecomposable tensor. Being written in components, this 
relation has the form 

C$X j U f , = YiV a , (6.3) 

where indices a, (3 (resp., i, j) correspond to C 2 (resp., 7i) and summation over repeated indices is 
implied. 



Suppose (6.2) holds; then the vector X is called a vacuum vector of the L-operator C. Multiplying 



(3.2) from the left by the covector V = (Va, — Vi), orthogonal to we get 

(VCU)X = . (6.4) 



Here VCU is an operator in TL with matrix elements V a C l !%Up. Conversely, suppose (3.4) holds. Then 



we have (3.2) with some vector Y which is uniquely determined by U, V and X. 

The relation ( |6.2|) (in the particular case TL = C 2 ) was the starting point for Baxter in his solution of 
the eight-vertex model [[hi} . In the papers on integrable lattice models of statistical physics this relation 
is called " pair propagation through a vertex" . In the context of quantum inverse scattering method 0] 



the equivalent condition (6.4) is more customary. It defines the local vacuum of the gauge-transformed 
L-operator (this explains the terminology introduced above). In general form this relation appeared for 
the first time in p3| . 

The necessary and sufficient condition for the existence of vacuum vectors is 

det(V£U) = 0. (6.5) 
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Putting for simplicity XJ% = V2 = 1 and using ( |6.l| ), one may represent ( |6.5| ) in a more explicit form: 



det(CM + B- U X V X C - V X V) = 0, U 2 = V 2 = 1. 



(6.6) 



This equation defines an algebraic curve in C 2 which is called the vacuum curve of the L-operator C. 
So vacuum vectors are parametrized by points of the vacuum curve (i.e. by pairs (U%,Vi) satisfying 
(6.6)). In general position the space of vacuum vectors corresponding to each point of the curve is 
one-dimensional. 

Suppose now that Ti = C 2 and C satisfies the equation ( fL"q ) with some matrix R. In this case the 
vacuum curve has genus 1, i.e., it is an elliptic curve £q. It is parametrized by points z of one-dimensional 
complex torus with periods 1 and r. 

Fixing a suitable normalization (for example, putting second components of all the vectors equal to 
1), we may consider components of the vectors U(z), V(z), X(z), Y(z) as meromorphic functions on £q 
having at most 2 simple poles. With this normalization, the right hand side of (6.2) must be multiplied 
by a scalar meromorphic function h(z). It follows p3| from the Yang-Baxter equation that 



Y(z) = X(z + ^), V (z) = U(z-^), 



where r\ and rf are some constants. Therefore, we can write the basic relation (3.2) in the form |33j 



(6.7) 



C(X(z) 2) U{z)) = h{z)Y(z) <g> V(z) = h(z)X(z + |) ® - y 



(6.8) 



Let Dx (resp., Djj) be the divisor of poles of the meromorphic vector X(z) (resp., U(z)). Let M.(D) 
be the space of functions associated to an effective divisor D, i.e., functions having poles at points of D 
of order not higher than the multiplicity of the corresponding point in D. For divisors of degree 2 on 
elliptic curves this space is two-dimensional in general position, so dim M.[Dx) = dim Ai{Djj) = 2, and 
components of the vectors X and U form bases in these spaces. Further, the functions Xi{z)U a {z) 1 i,a — 



1,2, form a basis in the space Ai(Dx + Du). According to ( 3.S), the functions h(z)Xi(z + Q)U a (z— y ) 
form another basis in this space and the matrix C connects the two bases. The divisors of poles of the 
left and right hand sides of (6.S) must be equivalent, i.e., must be equal modulo periods of the lattice. 



Since under the shift by rj' /2 the divisor of poles of the function having 2 poles is shifted by rj', this 
means that 

t 1 / -t 1 = M + Nt, M,NeZ. (6.9) 
The vectors X(z), U(z) are double-periodic, hence there are 4 different cases: 



rj = -q + 2w a , 



(6.10) 



where uj = 0, u)x = 1/2, uj 2 = (r + l)/2, lu 3 = r/2 (cf. (|5.26|) ). 
Baxter's parametrization of the L-operator follows from eq 



06, 



Indeed, the equivalence class of 
the pole divisor of X(z) may differ from that of U(z) by only a shift on £q. By means of a "gauge" 



transformation one may represent C in the form (1.4). The value of this shift is then i dent ified with the 
spectral parameter of the L-operator. In this parametrization, it is natural to write ( |6.8|) explicitly in 
terms of ^-functions. 

To do this, it is convenient to use another normalization, specifically, the one in which the vectors 
are entire functions in z (in this case they are sections of certain line bundles on £q). 

Let us introduce the vector 



e(z) 



V 



(6.11) 



2U 



Its components form a basis in the space of ^-functions 9(z) of second order with monodromy propertties 
6(z + 1) = 0(z), 6(z + t) = exp(— 2mr — 4iriz)6(z). They have 2 zeros in the fundamental domain of the 
lattice with periods 1 , r. Then, putting 

X(z) = Q(z), (6.12) 
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1 



U(z) = U ± (z) = Q{z±^(u+ 7 l) 
(for each choice of the sign) , we may rewrite ( |6.2[ ) in the form , 



LW(«)e(z)®e(*±-(u + i)) = 2 fl ± 



l(»7|t) 



_ ( -) (:+ 2) 00 (z±~( U -|)±O, 



where 



5? = 1, 



- exp(zb27rzz + 7ri(u ~\ )), 



7b+l 



(u\r) 



6=0 



W||r) 



(6.13) 

(6.14) 

(6.15) 
(6.16) 



(see dl.4| ) and (1.7)). We note that L^ a \u) = a a L(u) (L{u) is given by ( |l.4| ) with S a = a a and the matrix 
product is performed in the auxiliary space) satisfies the " RLL = LLR" relation (1.5) with the same 
i?-matrix ( |l.7| ) for each a = 0, ... ,3. The scalar factor in the right hand side of (6.14) is determined 
from the condition that £(§) is proportional to the permutation operator in C 2 (8> C 2 . One may verify 
(3.14) directly using identities for ^-functions (see the Appendix). 

Remark. Given an elliptic curve, we can always choose the vector X(z) to be an even function, 
X{—z) = X(z). We introduce the even function X(z) ( |6~lT| ) fr om the very beginning. Then the equality 
corresponding to the minus sign in (S.14) follows from the similar equality with the plus (it is enough 
to change z — > —z). However, in Baxter's approach it is useful to deal with both equalities. 



Let us turn to the case of arbitrary spin. Consider an L-operator of the form (1.4) 
L(u)-- 



W (u)S + W 3 (u)S 3 
Wi(u)S 1 +iW 2 {u)S 2 



W 1 {u)S 1 -iW 2 {u)S 2 
Wo{u)S - W 3 (u)S 3 



where 



W a {u) 



'a+l 



(u\t) 



WIIt) 



(6.17) 



(6.18) 



and S a are generators of some algebra (at this stage the commutation relations (pTl|), (1.2) are not 
imposed). We are going to obtain explicit formulas for representations of this algebra from a simple 
generalization of (6.8). 

Before presenting the main result of this section we need some more preliminaries. 

Since we are interested in the general form of representa tions in terms of difference operators, in 
what follows we take for the quantum space of the L-operator ( 6.17 ) the space of meromorphic functions 
of one complex variable z. The generators S a act on elements of this space (i.e. functions X(z)): 



S a : X(z) .— » (S a X)(z) 



Consider the following generalization of ( |6.S| ) and ( |6.14 ) 



U ± {z) T {a a L{u))X{z) = gaiu^iz T ivT to a ) T X(z ±j), a = 0, ... ,3. 

Here z £ £ , U T — (Ui, U 2 ), £ is a parameter, g a {u) are scalar functions independent of z, and 

u + £77 s 



u ± {z) = e[z± 



(6.19) 



.20) 



(6.21) 



(for £ — 1/2 it coincides with ( 6.13 )). As before , L(u) acts on U^ 1 as_2 x 2-matrix, while each matrix 
element of L(u) acts on X(z) according to ( 6.19 ). We have written ( |6.20| ) in terms of the covector U T 
for the following reason: since the operator (6.19) acts on a function from the left, this is equivalent to 
the right action of the corresponding matrix (representing this operator in a fixed basis) on the covector 
formed by components of the function with respect to the basis. As it is clear from what follows, at 
I — 1/2 ( 6.2C| ) coincides with the conjugated equality Q6.14 ). 
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It should be noted that at present time we can not suggest any explicit description of the vacuum curve 
of the L-operator ( |6.17 ). Moreover, we do not know any direct argument establishing the equivalence 
between ( |6.20| ) and the "intertwining" relation ( |1.5| ) for L(u) (taken together with the Yang-Baxter 
equation for R(u)). Theorem 6.1 (see below) states that Sklyanin's commutation relations (1.1, 1.2) for 
S a follow from ( |6.20 ), hence L(u) should satisfy (1.5). Let us stress once again that our arguments are 
in a sense inverse to original Sklyanin's approach (see also p4j , [ pq] , where some formulas for vacuum 
vectors in the higher spin XYZ model were obtained). Our starting point is the relation ( 6.20| ), where no 
conditions on S a are implied. It turns out that the Sklyanin algebra for S a (together with its functional 
realization) follows from (3.20). 

The main result of this section is the following 



Theorem 6.1 Let L(u) be given by (6.17), where S a are some operators in the space of meromorphic 
functions X{z) of one complex variable z. Suppose the relation (6.2C ) holds for a = 0, i.e., 



t/ ± (z) T L( W )X(z) = goiup+iz T ir,) T X{z ± |), 



(6.22) 



where U ± (z) is defined in ( 6.21 ), I is a parameter and go(u) is a scalar function independent of z. Then 
S a are difference operators of the following form: 



(S a X)(z) = X 



(»)'-'fl a+1 (f|T) 

0i(2z\t) 



8 a+1 (2z~£n\T)X(z+^) 



9 a+1 (-2z-er 1 \r)X(z-j) 



(6.23) 



where A is an arbitrary constant. Conversely, if S a are defined by (6.25), then (6.2i) holds, and 

g (u) = 2\6i(u + £r)\T). 



(6.24) 



Remark, a) Using transformation properties of the vector Q(z) under shifts by half-periods oj a it is 
readily se en th at the relations ( 6.20| ) for a ^ follow from ( 3.22| ); b) For even functions X(z) the two 
relations ( 6.22| ) are equivalent. 

Proof. The relations ( |6.22 ) form a system of 4 linear equations for 4 functions (S a X)(z) entering the 
left hand side. More explicitly, we have 



Uf{z)(L lx {u)X){z) + U 2 ± (z)(L 21 (u)X)(z) = g (u)U±(z T £ V )X(z ± i) 



(6.25) 



Ut{z)(L 12 {u)X)(z) + U±{z){L 22 {u)X){z) = go(u)U?(z T lrj)X{z ± j) 



(6.26) 



where (L a p(u)X)(z) are expressed throu gh (S a X)(z) according to ( |S.17 ). Fix go(u) to be given by 
( 3.24 ). Solving this system, we get ( 6.23|) (all necessary ide ntitie s for ^ -functions are presented in the 
Appendix to this section). Therefore, (6.23) is equivalent to (6.22, 6.24). ■ 



Corollary 6.1 Suppose for some L(u) of the form (6.17) the equation (6.2i) holds. Then L(u) satisfies 
the "intertwining" relation (1.5) with the R-matrix (l. r i). 

The proof follows from the identification of ( |6.23 ) with formulas ( |l.9| ) for representations of the Sklyanin 
algebra by putting 2z = x, X{x/2) = f(x). The constant A is not essential since the commutations 
relations (1.1, 1.2) are homogeneous. 



Remark. From the technical point of view, the derivation of formulas (1.9) by solving the system 



(S.22) is much simpler than the direct verification of the commutation relations. The amount of com- 
putations in the former case is comparable with that in the latter one if we identify the coefficients only 
in front of f(x ± 2rj). 



Let us consider now the equality (6.20) for a — 1, 2, 3. 



Lemma 6.1 Let us define the transformation of the generators Sb 1 
relation 

o a L(u + uj a ) = h a {u) V M y,^) g, ab , 



y a (S b ), 0,6 = 0, 



, 3, by the 
(6.27) 
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where h$(u) — h\{u) — 1, h,2(u) — —ih^iu) — exp(— ^-p — inu). Then y a , a = 0, ... ,3, is an automor- 
phism of the algebra generated by Sb- 



The proof follows from the fact that a a L(u + uj a ), a = 0, ... ,3 satisfies ( fL5| ) with i?-matrix (1.7). (This 
is because the matrices a a are c- number solutions of ( |l ■£) ).) From (6.27) we find the explicit form of y a 
(here O (§ ) = a (f|r)) : 



y\ : (So, Si, S2, S3) h 
3^2 : (5*0) Si, S2, S3) 
■ {So, Si,S2, S3) 



Ml), 



m) 



Sim °' 



^3(|) -no. 

~^3, n lr , \ ^2 



1(1: 



^(f) 



Ml) 
.M§) 

U(§) 



-S- 



Q V 2 / Q 



3(1) " 



Mf) q ^2 

9i(|) 

>->l 



L(i) 



So 



(6.28) 

(6.29) 
(6.30) 



and .X) is the identity transformation. These automorphisms were considered by Sklyanin in J8|. 

Shifting u — > u+cj a in ( 3.20 ) and applying these automorphisms to ( 6.23| ), we get for each 6 = 0, ... ,3 
the following representations: 



9i{2z) 



lz-£r]-u b )X{z + ^)-e a+l (-2z-£r)-w b )X(z ■ 



(6.31) 



We remark that (3.31) can be obtained from ( |6.23 ) by the formal change t — > t + LOb/rj provided £ is 
considered to be an arbitrary complex parameter. However, we will see soon that the operators (6.31) 
restricted to invariant subspaces yield non-equivalent finite-dimensional representations of the Sklyanin 
algebra. 

Suppose £ £ then one may identify X(z) with a section of some linear bundle on the initial 

elliptic curve £ . Repeating the arguments presented after eq. (|6.8|), we conclude that degree of this 
bundle equals At. From the Riemann-Roch theorem for elliptic curves it follows that in general position 
the space of holomorphic sections of this bundle is 4£-dimensional. It is convenient to identify these 
sections with ^-functions of order At. Consider the space T^t of even 9- functions of order At, i.e., the 
space of entire functions F(z), zeC, such that F{— z) — F(z) and 



F(z + l) = F(z), 

F(z + t) = cxp(-Atnir - Stiriz)F(z) . 



(6.32) 



In complete analogy with the paper B , one can see that the space is invariant under action of the 
operators (6.31) for 6 = and 6=1, while for 6 = 2 and 6 = 3 the invariant space is exp (—Tri z 2 /r])T 4 ^. 
It is known that dim 7^ = 2t+ 1 provided t G ^Z + . Restricting the difference operators (6.31) to these 
invariant subspaces, we get 4 series of finite-dimensional representations. 

Generally speaking, these representations are mutually non- equivalent. This follows from analysis of 
values of the central elements. The Sklyanin algebra has two independent central elements: 



(6.33) 



a=0 



(the constants Ji are defined by (1.3) and (1.10)). Their values for the representations (6.31) at 6 = 
0, . . . , 3 are: 

Ko=Ael(i£ + ^)r 1 + u b \^j, (6.34) 

K 2 = 46>i {(£ + 1)77 + w 6 |r) 6 t {tr, + oj b \r) . (6.35) 
The arguments given above lead to the following statement. 



Theorem 6.2 The Sklyanin algebra (1.1, l.i ) has 4 different series of finite- dimensional representations 
indexed by b — 0, 1, 2, 3. Representations of each series are indexed by the discrete parameter t £ 
("spin"). They are obtained by restriction of the operators (6.31) to the invariant (fit + I) -dimensional 
functional subspaces for 6 = 0, 1 and exp(— niz 2 /i~i)T 4 ^, for 6 = 2,3. For general values of parameters 
these representations are mutually non- equivalent. 
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Let us identify these representations with those obtained by Sklyanin in m. For 6 — and 6 = 3 
we get series a) and c) respectively. These representations are self-adjoint with respect to the real form 
of the algebra studied in The other two series (corresponding to b = 1 and b = 2) in general are 
not self-adjoint. Consider first the case 6=1. For rational values of i], i] = p/q, and special values 
of £, £ = (q — l)/2 mod q, these representations are self-adjoint and are equivalent to some subset of 
representations of series b) f\ To the best of our knowledge, the series corresponding to 6 = 2 was never 
mentioned in the literature (though, in a sense, it is implicitly contained in Sklyanin's paper). Another 
outcome of our approach is the natural correspondence between the different series of representations 
and points of order 2 on the elliptic curve. 

It is natural to surmise that representations of the last two series become self-adjoint with respect 
to other real forms of the algebra. A real form is defined by an anti- involution (^-operation) on the 
algebra. It should be noted that classification of non-equivalent real forms of the Sklyanin algebra and 
its generalizations is an interesting open problem. 



Concluding this section, we would like to remark that the variable z in ( 3.20 ) and ( 5.31 ) may be 
identified with the statistical variable (" height" ) in IRF-type models [|(| (after a suitable discretization 
of the former). This is readily seen from the well known vertex-IRF correspondence recalling that the 
vertex-IRF transformation is performed by means of the vacuum vectors (for the explicit form of this 
transformation in the case of higher spin models see |]35| ). 

Finally, it seems to be instructive to carry out a detailed analysis of the trigonometric and rational 
limits of the constructions presented in this section. Some particular related problems have been already 
discussed in the literature. In the recent paper [[37| , vacuum vectors for the higher spin XXZ-type 
quantum spin chains are constructed. Vacuum curves of trigonometric L-operators have been described 
in p8f . In the simplest case they are collections of rational curves intersecting at 2 points. Trigonometric 
degenerations of the Sklyanin algebra that are in a sense " intermediate" between the initial algebra and 
the standard quantum deformation of gli are studied in [^|) . 

Appendix to Section 6 

We use the following definition of the ^-functions: 

9x(z\t) = ex P + h 2 + 2m{z + h(k + h) , (6.36) 

fcez ^ ' 

Qi{z\t) = J2 ex P (*iT{k + i) 2 + 2mz{k + ~) j , (6.37) 
fcez V - / 

6s{z\t) = exp (irir k 2 + 2nizk) , (6.38) 
fcez 

8 4 {z\t) = ^ exp (-Kirk 2 + 2m{z + hk) . (6.39) 
fcez ^ ' 

For reader's convenience we recall here the definition of the er-function used in Sects. 2-4: 

^^=m^\^r^ Y (6 - 40) 

Here is the list of identities used in the computations. 
The first group of identities (addition theorems): 

6Mt)6Mt) = 0i(x + y^e^x - y\2r) + 6^ + y^O^x - y\2r), (6.41) 



3 The whole family of representations of series b) found in |8| has 3 continuous parameters. These representations are 
t ari d exist only if r\ = p/q. In this case all of them have dimension q. They are obtained by restriction of the 



self-adjoir 

operator (6.3f) to a finite discrete uniform lattice. 
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9 4 (x\T)e 4 {y\T) = 9 3 (x + y\2 T )6 3 (x - y\2r) - 9 2 (x + y\2r)9 2 {x - y\2r), (6.42) 

OMt)OMt) = 3 (x + V \2t)9 3 (x - y\2r) + 6 2 {x + y\2r)9 2 {x - y\2r), (6.43) 

6 2 (x\T)6 2 (y\T) = 9 3 {x + V \2t)6 2 {x - y\2r) + 6 2 (x + y\2r)9 3 {x - y\2r), (6.44) 

9 1 (x\T)9 1 (y\r) = 9 3 (x + y\2r)e 2 (x - y\2r) - 9 2 (x + y\2r)0 3 (x - y\2r), (6.45) 

Here are simple consequences of them convenient in the computations: 

9 4 {x\T)e 3 (y\T) + 9 i (y\T)9 3 (x\T) = 29 4 {x + y\2r)9 4 (x - y\2r), (6.46) 

OMtWMt) - 9 4 { V \t)9 3 {x\t) = 29 x (x + y\2r)9 x (x - y\2r), (6.47) 

6Mr)0 3 {y\T) + d 4 (y\r)9 4 (x\T) = 26 3 (x + y\2r)9 3 (x - y\2r), (6.48) 

9 3 {x\T)9 3 {y\T) - e 4 (y\T)9 4 (x\T) = 28 2 (x + y\2r)9 2 (x - y\2r). (6.49) 

The second group of identities: 

29M^)e A {y\2T)^e l {^\T)9 2 {^\T) + e 2 {^\r)9 1 (^\T), (6.50) 

29 3 (x\2r)9 2 (y\2r) = l( £±i| r )0 l( £Z^| T ) +0 2 (£+»| T )0 2 (£Z»| r)) ( 6 .51) 

20 3 (x\2r)e 3 (y\2T) = e 3 (^\T)e 3 (^\ T ) + e 4 {^\T)e4^\T), (6.52) 

2^(^12^(1,12^ =e 3 (^|r)6»3(^|r)-(94(^|r)fl4(^|r). (6.53) 
Particular cases of them: 

29 1 {z\r)e 4 (z\r) = eM T -)9 1 {z\\), (6-54) 

29 2 {z\t)9 3 {z\t) = 9M T -)9 2 {z\ T -). (6.55) 

Two more identities: 

e 1 (z\^)9 2 (z\^)=9 i (0\r)e 1 (2z\r), (6.56) 

8M\)0M\) = 4 (0\t)9 4 (2z\t). (6.57) 



7 Concluding remarks 

This work elaborates upon the following three subjects: 

I) Dynamics of poles for elliptic solutions to the 2D non-abelian Toda chain; 

II) Difference analogs of Lame operators; 

III) Representations of the Sklyanin algebra in terms of difference operators. 
Let us outline the results: 

- The poles move according to equations of motion for spin generalizations of the Ruijsenaars- 
Schneider model; the action-angle variables for the latter are constructed in terms of some algebraic- 
geometrical data; 

- One of the generators of the Sklyanin algebra, represented as a difference operator with elliptic 
coefficients, has the "finite-gap" property that is a motivation for the analogy with Lame operators; 

- Starting from the notion of vacuum vectors of an L-operator, a general simple scheme for construct- 
ing functional realizations of the Sklyanin algebra is suggested. 

Here we would like to explain why the three themes are to be intimately connected. 
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To each problem I) - III), a distinguished class of algebraic curves has been associated. In case I), 
these are spectral curves T for the L-operators of the Ruijsenaars-Schneider-type models; in II), we deal 
with the spectral curve T' for the difference Lame operator So (a generator of the Sklyanin algebra); 
in III), the representations are defined on sections of certain line bundles on a vacuum curve £ of the 
elliptic higher spin L-operator ( |1.4| ) (though it is implicit in Section 6). It has been shown that T and 



r' are ramified coverings of the initial elliptic curve. The characteristic property (6.20) of the vacuum 
vectors suggests that the same should be true for £ , i.e., £ is a ramified covering of the initial elliptic 
curve £q (the vacuum curve of the spin-1/2 L-operator). 

The connection between I) and II) is similar to the relation between elliptic solutions of KP and KdV 
equations. Specifically, the elliptic solutions of the abelian 2D Toda chain, which are stationary with 
respect to the time flow t + + f_ correspond to isospectral deformations of the difference Lame operator 
Sq (considered as a Lax operator for ID Toda chain). In other words, the hyperelliptic curves T' form a 
specific subclass of the curves T. A similar reduction in the non-abelian case yields spin generalizations 
of difference Lame operators. Their properties and a possible relation to the Sklyanin-type quadratic 
algebras are to be figured out. 

Apart from the apparent result that the construction of Sect. 6 provides a natural source of difference 
Lame-like operators, we expect a more deep connection between II) and III). Specifically, the spectral 
curves V are expected to be very close to the vacuum curves £ . Conjecturally, they may even coincide, 
at least in some particular cases. At the moment we can not present any more arguments and leave this 
as a further problem. 



At last, we would like to note an intriguing similarity between the basic ansatz ([2.26 ) for a double- 



Bloch solution of the generating linear problem and the functional Bethe ansatz ]39j. Indeed, in the 
latter case wave functions are sought in the form of an "elliptic polynomial" T[ °~( z ~ z j), where the roots 
Zj are subject to Bethe equations. Similarly, in the former case we deal with a ratio of two "elliptic 



polynomials" (cf. (5.3C )). However, this function is parametrized by residues at the poles rather than 
zeros of the numerator. This may indicate a non-trivial interplay between Calogero-Moser-type models 
(and more general Hithin's systems) and quantum integrable models solved by means of Bethe ansatz. 
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